Chiral violations in domain-wall QCD 
from one-loop perturbation theory at finite A^^ 



Stefano CapitaiiQ 
Institut fiir Physik, FB Theoretische Physik 
Universitat Graz, A-8010 Graz, Austria 

We present perturbative calculations made with domain-wall fermions which possess a finite num- 
ber of points Ns in the extra fifth dimension. We have derived the required propagator functions, 
investigated the one-loop properties of quark amplitudes at finite Ns and evaluated three quantities 
that can provide insights on chirality-breaking effects from the perturbative side. 

First we have computed the residual mass for various choices of Ns and of the domain- wall height 
M. We have found that this radiatively induced mass approaches zero reasonably fast with the extent 
of the fifth dimension, depending on AI and on a lesser extent on the coupling go. We have also 
computed the differences of the renormalization constants of the vector and axial-vector currents 
and of the scalar and pseudoscalar densities. Finally we have calculated the chirally-forbidden 
mixing (which at finite Ns is suppressed only partially) of an operator which describes the lowest 
moment of the g2 structure function. In general we see that at M = 1.8, where simulations are 
usually performed, values of Ns = 20 or larger would be desirable in order for chiral violations to 
be negligible. 

The quantities that we have studied turn out to lose gauge invariance when Ns is not infinite. 
We have also found that anomalous dimensions of operators at finite Ns generally depend on Ns 
and M. In particular, the vector and axial- vector currents have in general a nonzero anomalous 
dimension at finite Ns- 

PACS numbers: 12.38.Gc,11.30.Rd,11.30.Qc,11.10.Gh 



I. INTRODUCTION 

Domain-wall fermions Q, |^ constitute one of the 
known solutions of the Ginsparg- Wilson relation Q and 
are hence invariant with respect to chiral symmetry 
transformations even away from the continuum limit, for 
nonvanishing values of the lattice spacing a 0, while 
avoiding at the same time unpleasant effects like doublers 
and nonanaliticities. The massless chiral mode generated 
by the domain wall survives 1-loop renormalization and 
the theory has been proven to be renormalizable at this 
order Q . The possibility arises of the construction 
of chiral gauge theories at finite a Moreover, with 
this kind of fermions the leading discretization errors are 
reduced to O(a^) when chiral symmetry is exact. 

A certain amount of chiral symmetry breaking arises 
however in Monte Carlo simulations of domain-wall 
fermions, because they must be performed using lat- 
tices which have a finite number of points, Ng, in the 
fifth dimension. It is only in the theoretical limit in 
which the extension of the fifth dimension becomes infi- 
nite that the chiral modes (which are exponentially con- 
fined on the two opposite walls) can fully decouple from 
each other, yielding an exact chiral symmetry. The chi- 
ral modes acquire some mass if the distance between 
the two walls is not infinite, and to study the mag- 
nitude of these chirality-violating effects is one of the 
main objectives of the present work. With the com- 
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puter speeds presently available it is unfortunately not 
yet possible to obtain significant physics from simula- 
tions performed at a large Ng, where these chiral vio- 
lations would be numerically negligible. After the first 
pioneeringMonte Carlo implementations of domain-wall 
fermions | g , UOl and su bseq uent advances reported in 
[H III m IM El and [II 113, the most recent sim- 
ulations, which have considered phenomenological quan- 
tities as diverse as weak matrix elements, structure func- 
tions and heavy-light meson spectroscopy, have been 
mostly performed using lattices with only Ns — 12 or 
16 (for a selection of the latest results see for example 

IS lii IS il il 111 111 ii 111 111 111 1^ 

The experiences gained in these recent investigations 
seem to indicate that for extents of the extra fifth di- 
mension as small as Ns = 16 the chirality-breaking ef- 
fects Hi mm mm Hi are stm under control. The 
residual mass rUres for typical lattice spacings of about 
= 1.5 — 2 GeV generally turns out to be of 0(10^^) 
or O(10~^), depending on whether simulations are car- 
ried out in quenched or full QCD, and on which type 
of gauge action is used 

Ellllllll- In particular, when 
using quenched QCD instead of the full theory, or renor- 
malization group improved gauge actions (like Iwasaki 
and DBW2) instead of the simple plaquette action, the 
residual mass becomes substantially reduced. This is re- 
lated to the fact that, for a fixed choice of a, the values 
of (3 = are larger in those cases, and the gauge 

fields smoother. The residual mass is then in general 
not very small (expecially in the full QCD case), and it 
can at times become comparable to the input light sea 
quark masses. However, if the exponential suppression of 
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chirality-breaking effects takes place rather quickly with 
the length of the fifth dimension, then increasing Ng a 
little further could already be sufficient to obtain at last 
almost negligible chirality-breaking effects. 

The extent to which chiral symmetry is broken in lat- 
tices with a small Ng is thus one of the most important 
issues which need to be understood in present domain- 
wall simulations. To the extent that one- loop calculations 
can provide clues to the true behavior of the truncated 
domain- wall theory at finite Ns , we believe that it is use- 
ful and interesting to investigate chirality-breaking effects 
also from the point of view of perturbation theory, com- 
plementing nonperturbative investigations of such effects. 
Towards this end we present here the results of some se- 
lected one-loop calculations which we have carried out 
using the Feynman rules which exactly correspond to the 
theory at finite Ns . This is at variance with past domain- 
wall perturbative calculations where, in place of the ex- 
act quark propagators, their asymptotic expressions for 
large Ng were instead used [33, 111, |33, 113 • The purpose 
of this work is to calculate with the exact Feynman rules 
the deviations from the = oo results in the case where 
Ns is limited to small values of 0(10), that is for situ- 
ations which roughly correspond to present simulations. 
We use the plaquette gauge action, and we can obtain 
some estimates of the amount of perturbative chiral vio- 
lations by focusing on quantities like the additive renor- 
malization to the quark mass (i.e., the residual mass) 
and the deviations of some amplitudes from their values 
at A'^s = oo, including a chirally-forbidden mixing. Since 
the cost of domain- wall simulations grows approximately 
linearly with Ns , it is of some importance to understand 
how small Ng can be kept without occurring in large val- 
ues of the residual mass. A thorough exploration of large 
regions in the two-dimensional space spanned by Ng and 
the domain- wall height AI would be quite expensive when 
using Monte Carlo simulations, and perturbation theory 
remains then often the more practical way for gathering 
hints of what is happening in this space. The study of the 
dependence on M and Ns of various indicators of chiral 
violations is thus one of the main aims of the work. 

Significant perturbative insights at finite Ns be- 
yond tree level have been provided some years ago by 
Kikukawa, Neuberger and Yamada who diagonal- 
ized the mass matrix in the truncated overlap and de- 
rived 1-loop equations under certain assumptions on the 
gauge fields. We improve here on this by providing nu- 
merical results for the residual mass. What we present 
is a complete 1-loop calculation of such radiative effects, 
with the aim of seeing how mres, and other quantities 
which can act as indicators of chiral symmetry breaking, 
behave when Ns and M change. 

Other methods have also been used to gain insights on 
these effects. Recently Christ ^3 investigated the 
residual mass analytically by looking at the eigenfunc- 
tions of the five-dimensional transfer matrix [3| . Building 
on the understanding of the localization properties of the 
domain- wall modes, characterized by their mobility edge 



Ac |4a,lM|45| (for recent investigations see |4g . l47ll4S|| '). 

the leading effects were estimated as 



exp(-AcjVs) 
Ns 



RfpiiO)^, (1) 



where p is the density per unit spacetime volume of the 
eigenvalues of the logarithm of the transfer matrix, and I 
and e stand for localized and extended modes with aver- 
age size R respectively. The energy threshold from local- 
ized to extended modes is given by the mobility edge Ac, 
which is then also responsible for the speed with which 
the chiral violations decay with Ng . 

This article is organized as follows. In Sect, ^we in- 
troduce all propagator functions which are necessary for 
the perturbative calculations at finite Ns, and in Sect. 
mil we analyze the one-loop renormalization of quark am- 
plitudes in this theory. Since in this article we restrict 
ourselves to the calculation of finite diagrams, the treat- 
ment of divergences in the theory at finite Ng is left for a 
future work. We have however computed the coefficients 
of the divergent terms in a few cases, and we have seen 
that they are in general not equal to their continuum val- 
ues. In fact, they depend on Ng and on the height of the 
domain wall, M, as we show in Sect. IIVI where we also 
briefiy discuss the implications of this finding. In Sect. 
Ivl we then present the computation of the residual mass 
and encounter another feature of calculations at finite 
Ng, namely that gauge invariance is lost. In Sect. IVII we 
show the results at finite Ng for the difference between the 
vector and axial-vector renormalization constants, which 
should be zero at infinite Ng. In Sect. IVIII we finally 
present a power-divergent mixing due to the breaking of 
chiral symmetry, for an operator which describes polar- 
ized parton distributions, and in Sect. IVIIII we discuss 
what happens near the borders of allowed values of M, 
M and M — + 2, before making in Sect. IIXI some 
concluding remarks. 



II. PERTURBATION THEORY 

We work with the standard formulation of domain-wall 
fermions of Shamir 



S. 



DW 



X s—1 L M 

-V's(a;)(7M + r)Ul{x - fi)ips{x - //)) 
+ (i^six)P+^s+i{x)+7jg{x)P-ijg-i{x)') 



+ {M - l + 4r)tljs{x)ips{x) 



(2) 



X 

where we put r = — 1 , that is the Wilson term is added to 
the action with minus the conventional sign. The height 
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of the domain wall, or Dirac mass, M, at tree level sat- 
isfies < M < 2, so that the correct pattern of chiral 
modes (with no doublers) is attained when ^ oo: in 
fact for M < there is no chiral mode, while for M > 2 
there are four of them (and even more when M is further 
increased). The chiral projectors are P± = (1 ± 75)/2. 
Here and in most of the paper we put a = 1, but in 
some contexts, like when discussing the residual mass, 
the lattice spa cing will be explicitly shown. 

We refer to |38l Is^ l40l | for the Feynman rules which 
derive from this domain- wall action in the limit — > oo. 
In this article we always work at finite N^, and the ex- 
pressions that we have to use for the quark propagators 
are then different. They were partially derived in |37| |. 
and here we compute the remaining functions and pro- 
vide the complete set of required propagators. 

In (four-dimensional) momentum space the domain- 
wall Dirac operator has the form 

Dst{p) = (5s,t ^i7pSinp^ (3) 

+ {W+{p) + mM+)P+ + {W-,{p) + mM-,)P^, 
where the mass matrix is given by 



In more explicit form, we have 



^7t = ^s,l6t,Ns, 



and 



W{p) = l-M-2r^sin 



;„2 PA 
2 ■ 



(4) 
(5) 
(6) 



(7) 



VF+(p) 



W-{p) 



/-W{p) 



(~W{p) 

1 



-W{p) 



-W{p) 



-W{p) I 



-W{p)J 



, (8) 

, (9) 

(10) 
(11) 



In practical terms this theory looks like having several 
flavors of lattice Dirac fermions, which are mixed in a 
very special way so that a large mass hierarchy is gener- 
ated. At the end this theory indeed contains one chiral 
mode which is nearly massless together with iVg — 1 heavy 
fermions. 

In this work we only consider the case in which no ex- 
plicit mass term appears in the Lagrangian (m = 0). The 
tree-level 5-dimensional quark propagator is then given 
by 



(^,(-p)^,(p)) = E [( - ^in^'A^ + ^suip) ) G^tip) P+ + {- sinPM ^s,u + W+{p) ) G^P) P- , (12) 

U 

where the expressions of the functions G^{p) and G^{p) are 



Gftip) 



Mp) 

F{p) 



(1 - VF(p)e^"(f))(e-2^="(f) - l)e("+*)"(P) + 2Wip) sinh(a(p))(e(''-*)"(P) + e-^''-*)"^^)) 



-(1 - M^(p)e"(P))(l - e2JV="(p))e-(''+t)"(p) 



(13) 



G^tip) 



A{P) 

Hp) 



(e-2"(p) „^(p)e-"(p))(e-2A'="(p) - l)e(''+*)"(P) + 2W^(p) sinh(a(p))(e('^-*)"(P) -f- e-^'^-*)"'^)) 



+ (e2"(p) _ iy(p)e"(p))(l - (PNMp)y-~(s+t)o.(v) 



A(p)(e(^=-l^-*l)"(f) +e-(^=-l^-*l)"(f)). (14) 



In these expressions the quantity a{p) appears, which is defined by the positive solution of the equation P,|4 

I l + iy'(p)+EAsinVA 
cosh(a(p)) = ' (15) 
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and one uses the abbreviations 
1 



Aip) 
F{p) 



2W{p) sinh(a(p)) 2smh{Nsa{p))'' 



(l-W^(p)e"(P)) 



(16) 
(17) 



.e-^="(p)(l-W^(p)e-"(p)). 



When W is negative, a situation which arises only when 
1 < M < 2 if the momentum is small enough, the prop- 
agator is given by the above equations with the replace- 
ments 



W 



'\W\, 



(18) 
(19) 



which imply that also sinh a changes sign. 

The standard "physical" quark fields used both in 
Monte Carlo simulations and in perturbative calculations 
are given by 

g(x) = P+Vi(x)+P_V^^^(a;) (20) 
g(x) = ^^^(x)P++^i(.t)F_. (21) 

Strictly speaking these fields do not correspond exactly 
to the chiral modes, which should be instead eigenvectors 



of the mass matrix, like 



(22) 



where one calls 



Wo = W{0) = 1 - M. 



(23) 



The physical quark field q{x) is however more convenient 
to use than xoix), given that wq undergoes already at 
one loop a renormalization deriving from the additive 
correction to the domain-wall height M (see Eq. H52f) 
below). Thus one instead takes q{x) to represent the 
physical zero modes of the theory. Moreover, at finite 
Ns an additional issue about xoix) would arise, because 
this field is an eigenvector of the mass matrix only up to 
terms of order Ns e-^'"(o) [Ullilill. 

The computation of matrix elements involving phys- 
ical states and operators requires the introduction of 
additional propagators which connect the 4-dimensional 
physical quark fields with the 5-dimensional quark fields 
which appear in the Lagrangian. We have here derived 
the expressions of these propagators for the case of finite 
Ns, and they are given by 



{q{-p)i^s{p)) - P+{M-p)i's{p))+P-{i^NA-p)i^s{p)) 

( ilt^sinpf, -nMp) 2Ty(p) sinh(a(p)) 
V E{p) ^ S(p)(l-e^2JV,"(p)). 

X (^(e"^^^"''^"^^^ - e-2W.a(p)g(JV.-s)Q(p)^p^ ^g-(s-l)a(!3) _ g-2Ar,a(p)g(s^l)Q(j9)jp_j 

^ g-"(p) ( ('g-(s-l)a(p) _ g-2(JV,-l)Q(p) (5-1)q(p)\ p /24^ 

_^(^g-(JV.-s)Q(p) _ ^-2iN,-l)a{p} ^(N,^s)a{p)^^ p_y 

(M-pMp)) = {^s{-p)i^Ap))P- + {M~p)i^NMP+ 

= |^(e~(^=~^^"(^') - g-2^="(P)g(^=-'')"(P))P_ + (|g-(s-l)"(p) _ g-2Af=a(p)g(s-l)a(p)^p^^ 

^ f il^sinpf, ^-N,a{p) 2T/F(p)sinh(Q;(p)) \ 
V E{p) £;(p)(l-e-2A^="(p))/ 

^ g-"(p) ( /'g-(s-l)a(p) _ -2(jV,-l)Q(p)g(6-l)Q(p)\p_ /25^ 

_|.(^g-(JV.-s)Q(p) _ g-2(Ar,-l)Q(p)g(A',-s)a(p)-jp^^ ^ 

where we have defined 

E{p) = 1 - W^(p)e"(P) - g-2Ar,a(p)^^ _ ^(p)g-"(p)). (26) 

We have also numerically checked the above expressions and the correctness of their implementation in our computer 
codes by verifying the validity for each s (for various choices of Ng) of the identities 

GUp)=G^i{p) = --^ (e-("-i)"(f) -e-2^="(P)e("-i)"(f)), (27) 
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^{p) = G«v ip) = (e-(^-^)"(p) - ^-2NMp)^iN.-sMp)\ (28) 

t t 

I p-A'=a(p) 2VF(p) Sinh(a(p)) ( (s^l)a{p) _ -2N^,aip) (s-l)a{p)\ 

£;(p)(l-e-2^="(p)) ^ ^ 

J2wrdp)G^siP)-T.^Jtip)GfAp) = - ^ _ e-"(rt _ ,-2(A^.-i).(p)^(.-i).(rt) (30) 

t t 

, _-Ar,a(p) 2VF(p) sinh(Q(p)) . (jv,_^)a(p) _ -2Ar,a(p) (Ar.-s)a(p)\ 

^ E{p) {1 - e-^^Mp)) ^' 

which relate the propagators {qi—p)4'sip)) {'4's{—p)q{p)) to {il;s(—p)ipt{p)) ■ 

The calculation of perturbative amplitudes also requires the knowledge of the expressions of these new propagators 
for small momentum. In this limit we obtain 

1 - " p'' +Wq (1 - Wq) ^ ^ 



where is defined in Eq. Since = e 

is easy to see that the terms which are proportional to 



^-N^aio) j-apidly approach zero when Ng becomes 
large. In the derivation of the above formulae we have 
used the useful small momentum expansions 



1 - iy(p)e"(P) 
1 - VF(p)e-"(f) 



P 



(33) 



= 1 



^0 1 _ ,,,2 P (34) 



and 



2P^(p)sinh(a(p)) ^ <°(l-^§) 



E{p) 



2 I 2Afs 

p^ +Wq 



. ^^^^ 

and we have dropped all terms of order p^Wp " and higher, 
which are much smaller than either of the factors p^ or 



alone, and are not relevant when p — > 0. It is easy 
to check that all the propagators that we have derived 
in this Section reduce to the expressions used in the cal- 
culations of Refs. [3^ ^ when A''^ is large. Notice 
also that the function E(p) introduced here tends in this 
approximation to the function F{p) as defined in those 
articles. The above expressions for the small momentum 
propagators can also be used for M > 1 without further 
modifications (as we have also numerically checked). 
Finally, we also need the function that describes the 



propagation of the physical fields alone. This is given by 
{q{-p)qip)) = ^{^l,8inp,il-e-^^^'^(P^) (36) 



-nMp) .2W{p)smhia{p)) 



which in the limit of small momentum becomes 



{q{-p)q{p))c 



(i-^g) i^,^?;!!""^ - (37) 



pi 



2N, 




It is interesting to see that it is also possible to calculate 
{q{^p)q{p)) in alternative way directly from the 5- 
dimensional propagator of Eq. I|12|l : 

(qq) = P+{Mn^P+ + P+{Mi)P- 

+P_ (V'w>^jP+ + P_ {'4'nJji)P- 
1 



sinp^ (g^^j^^ P+ + Gfi P- 



ly G^V^P+ + G^^lP_ 



(38) 



The fact that in this way we obtain again the result of 
Eq. l|H7|) provides a good check of the above formulae. 

Domain-wall fermions present thus at finite some 
new peculiar features. Although the theory which we 
have started from is described by a Lagrangian of mass- 
less quarks, the propagator of the physical quark field. 
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{q{—p)q{p))c, acquires when Ng is kept finite a nonvan- 
ishing mass term, which at tree level is given by 



,(0) = -,„^= 



(l-w^) = -{l-MY^" M{2-M), (39) 

as can be seen from the general expression of a fermion 
propagator of mass fi for small momentum in Euclidean 
space: 



(40) 



We readily see that this tree-level residual mass vanishes 
when Ns = oo. We will only consider even values of 
Ns , in which case the fermion determinant can be proven 
to be positive (so that the square root of the two-flavor 
theory is well defined and an odd number of dynamical 

flavors can be simulated) . Then mlel is always a negative 
quantity. With our calculation we have thus reproduced, 
up to a sign, the result for found in '5', '32', '33', "a?, "iTl . 
This was derived by considering the quadratic operator 

D, which could perhaps explain the sign discrepancy. 

We will see that when the one-loop corrections are 
taken into account, the residual mass changes sign and 
becomes positive. 

If we had used the chiral mode xo, Eq. (|22|) . the prop- 
agator for small momentum would have been given by 



Xo deviates from what would be the exact chiral mode 



for finite Ng by terms of order Ng 



The real physical 



propagator would then also present further corrections 
terms, and it could be that the expression in Eq. I|41|) 
would get simpler. 

To complete the setup of our calculations, we also recall 
that we use the plaquette action in a general covariant 
gauge, where the gluon propagator is given by 



G^.(fc) 



4E 



• 2 k„ 



(l-a) 



4 sin % sin ^ 
sin^ ^ 



4E. 



(42) 

so that a = 1 and a = correspond to the Feynman and 
Landau gauges respectively. The QCD vertices that we 
need are the usual ones, and (apart from color factors) 
they have the form 



-90 «7m ' 



(P) 



1 



' 2 

,^ 
2 



cos 



2 
2 



(43) 
(44) 



for the interaction of the quark current with one gluon 
and two gluons respectively, where p in this case stands 
for the sum of the incoming and outgoing quark mo- 
menta. 



((1 



2\ ^ 



XoXo)c 



l-wl 



2N, 



..Ns 



(1 



(l-z«o')(l-2u;^^=) 



-Ng 



(1 



2Ns 



(41) 



which compared to {qq)c has correction terms of higher 
order in Wr, " in the numerator. We remind however that 



III. RENORMALIZATION AT FINITE iV^ 



Let us now investigate the properties of the renormal- 
ization of the self-energy of a massless quark at finite A^^- 
We first notice that the one-loop propagator of the phys- 
ical field can be written, according to the general form of 
the external legs (g(-p)V's(p))c and {ips{-p)q{p))c, as 



{q{~p)q{p))i loop = 



1 



■Ji 2\ ' N,/^ 2\ '^I^P^ • / TVs /I 2^ 

tp — Wq [1 — Wq) tp — Wq [1 — Wq) "" ""^ 



=(1 - wl) 



^Ji~w^^{l~wl)~{l~wl)^,{p)' 



(45) 



where 
S,(p) 



s,t , 2Ns 

1 - Wn 



-Wo 

■^st{p) 

1 

I — w, 



N, 



27V, ~(Ns-s)\ r, 



2Ns 



Wo 'Wq 



, s-1 2Ns -(s-l)\ 

[Wn -Wq "Wq^ ')P- 



2(Ns-l) -(s-l)\ 



w-^- -'w^- ->)P+ + {w^^-^-wl^''^-'W'''~'^)P- 



(46) 



)P^ + {wl-'-wl^^w-^'~'^)R 
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-wo[[Wq ~Wq' 'wq' ')P- 



/ 1 (Uq 



contributes to the order and the functions and G^j give contributions to Sgt (p) only. The general form of Eg (p) 
is 



1-^0 



So 



— +¥[c'^, logo P + Eij - (i^- =(1 - Wo)) Y— ^ S3 



(47) 



where from now on we call for brevity = (^q/IBtt^) Cp 
(with Cf = {N^ - l)/2Nc for the SUiN^) gauge group). 
We can observe that differs in a few aspects from the 
corresponding expression for infinite Ng. In fact, apart 
from the different coefficient of the logarithmic term, 
which now depends on Ns and AI (see Section ITvTl . and 
the slightly different coefhcient of E3, it also contains a 
totally new contribution proportional to 1/a, called Eq, 
which is a mass correction term. The one-loop radiatively 
induced mass is indeed given by 



a m. 



(1) 



w^(l 



-2 

9 So, 



(48) 



as can be easily seen when the one-loop correction to the 
quark propagator is cast in the same form as its tree-level 
expression; 



{qi-p)<iip))i loop 



1 



¥ - (1 - wl) - (1 - wl) Sg(p) 



Wo 



(1) 



(49) 



This 0(a"^) critical mass, reminiscent of the analogous 
quantity for Wilson fermions, vanishes when the theory 
describes exact chiral fermions, that is at infinite N^, 
but is different from zero when computations are done 
at any finite Ng. This means that in the latter case the 
So contribution generates a finite additive renormaliza- 
tion to the quark mass, which can represent a measure 
of chirality-breaking effects. We associate this perturba- 
tive critical mass rrires with the residual mass which in 
Monte Carlo simulations is computed by looking at the 
explicit chiral symmetry breaking term in the axial Ward 
identities. 

The results above should not come as a surprise, since 
after all we are working here with a theory of Ng Wilson 
fermions, A^s — 1 of which are heavy states (not count- 
ing the doublers), and the necessity of a fine tuning to 
some critical mass comes out naturally. As in the case of 
Wilson fermions, where the hopping parameter is renor- 
malized away from its tree level value 1/8, in order to 
obtain a massless pion in domain-wall simulations the 
quark mass must be tuned to a nonzero number, which 
at one loop is given by the critical mass given above. 
This defines the chiral limit when no explicit mass term 
appears in the Lagrangian. Of course higher loops and 
nonperturbative effects give further contributions to the 
shift of the critical mass. From a practical point of view. 



it is interesting to see how small this critical mass is, to- 
gether with its dependence on Ng (and M). In the free 
case the numerical results for rrires according to Eq. H39() 
are collected in Table (where they have already been 
multiplied for IStt^, so as to ease the comparisons with 
the one-loop results for the critical mass presented in the 
next Section). We expect that the critical mass substan- 
tially decreases when Ng becomes large, as can be verified 
at one loop from Tables IXIll to IXVl in the next Section. 

The other quantities appearing in the last line of 
Eq. igU) are 



2^2 = 1 + .9 ( log a p -f E] 



(50) 



which is the quark wave function renormalization factor, 
and 



Z,,, — 1 



2W0 _2 ^ _ -, , ^2 



■g S3 = 1 g Zn 



(51) 



which represents an additive renormalization to wo and 
hence to the domain- wall height M Is^, as can be in- 
ferred from 

(1 -wl)Z^ = l- {wo + f S3)' + 0{f). (52) 

There is indeed no chiral symmetry which can protect 
this mass, even at Ng ~ 00. The additive renormal- 
ization to M can be traced back, in the damping fac- 
tors {q{-p)ipg(j)))c and {ips{-p)q{p))c in Eq. to the 
terms which are proportional to i/. Notice how these 
terms are also proportional to Wo'{^ — Wg), which al- 
though being of a different order in a is required for the 
correct recasting of the one-loop propagators in the form 
of Eq. P^ . 

The renormalization of a composite operator 
q{x) O q{x) which is multiplicatively renormalizable 
can also be expressed in a simple way. Again, by looking 
at the general form of the propagators the one-loop 
matrix element of such an operator between "physical" 
quark states can be written as 



{{qOq)qq)i loop 



1 



1 



Ao{p)-0 



Wo 



(53) 



where Ao{p) contains the contribution of the damping 
factors. For a logarithmically divergent operator it takes 
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the form 

Ao{p)^f{-j^o'''''^logaV + Bo), (54) 

where the anomalous dimension turns out in general to 
be a function of Ng and M, even at lowest order. 



erator, Eq. (|54|l . the coefficients of the logarithmic diver- 
gences turn out to depend on and M. It is only when 
Ns — oo that they become equal to the ones calculated 
in the continuum. We give here some examples of this 
phenomenon by computing a few of these coefficients. 



IV. DIVERGENCES AT FINITE A^, 

The divergence of the half-circle diagram of the sclf- 
As we have anticipated in the previous Section, in the energy comes from the terms which are of first order in 
case of the self-energy, Eq. (1471) . and of a divergent op- p in the integral 



2i 



2Ns 



(55) 



k ■ p 



X (^6px - (1 - a) ■ ^''^^ 



{k-pY 



{k-pf 



Ip ■ J^k^ ■ IX ■ (Gf,(fc) P+ + Gfi(fc) P_) 



2N^ -{N^-t) 



)p„-i-K-i-^r=<(*-^y^.), 



where the small k expansions of the functions and are given by 



2Ns 



2N,-s-t 



1 — W, 



2N^ 



-{s-t) 



..2N^ 



G%{k) 



w, 



s+t-2 



1 — W, 



2Ns 



- W, 



-(s-t) 



+ o(fc')j , 
)+o(fc2)), 



(56) 
(57) 



The terms proportional to 



2Ns-s-t 



and w, 



s+t-2 




which are localized near the two walls and provide the leading 



approximation for large Ng to the continuum coefhcient, have already been given in [stI l38| . 

After doing the gamma algebra and carrying out the sums in the indices s and t in the fifth dimension, we get a 
compact analytic expression for the coefficient of the logarithmic term, as a function of and M: 



(N.,Af) _ o 



1 — K7i 



2Ns 




1 - TVs 



1 — Wi 



2N^ 




2(Af,-|-l) 



1 — w. 



2N^ 



l-Ng 



l-wl 



2 n 2Ns \ 
Wa(l -w„ =) 



(58) 



where = a is the value of the coefficient in the case of exact chiral symmetry. Numerical values of c^^""*^-* for 
various choices of Ng and M in Feynman gauge are shown in Table ITU 

For a generic bilinear q{x)T q{x), the divergence is obtained by computing the following integral from the vertex 
diagram: 



d^k (1 - 



2N^ -(Ws-s)^ 



Wn ' W, 



2{N,-l)^^-(s~l} 



(2^)4 (l-«;g^=)4 ^paEs* 



2N, 



P_ 



(59) 



Wn Wr, 



P_ 



/ s-1 2Ns r> 

[Wn - Wn "Wn' ' P I 



■U _ (0) 



Wq 



1-wl 



fc2 + {mitf 



xF- 



^(^pA-(l-a)^^^) 



(k-p) 



(0) 
TTlres 



k^ + (mit) 



(K 



-Wn "W' ']P+ + (Wn -Wn "w' ' ) P- 
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The coefficient of the logarithmic term turns out after many simphfications to be given by 



2(^.-1) -(^.-m 



4 ' = cr • 1 - A^. u^o m: 1 - ^« ^0 2jv: - 4^^o U - 2]v: ' (60) 

where is its continuum value, that is = — —3 + a, Cy — c^^ = —a and — 1 ~ a. The dependence on A^s 
and M is different from the one of the self-energy, and thus we find that the anomalous dimensions of all bilinears 
must also depend on these parameters. Numerical values of c'y"'^^ are reported in Table Hill for the Feynman gauge. 

A remarkable consequence of the above formulae is that the anomalous dimension of the vector and axial-vector 
currents does not vanish anymore. The contributions coming from the half-circle diagram of the self-energy and 
from the vertex diagram of the vector (or axial- vector) current do not compensate each other, because part of the 
subleading terms are different. In Feynman gauge one gets 

7f = 2 wr^^ -{l-N. ) I { 2 + ) - 2 - ) , (61) 




and the numerical values of this anomalous dimension 
are reported in Table IIVI These currents have then for 
any finite A^^ an anomalous dimension which is different 
from zero. Only in the Landau gauge it is equal to the 
case of exact chiral symmetry and thus vanishes, but this 
happens just because the coefficients of the logarithms, 
being proportional to their values at A^s = oo, in this 
gauge vanish separately for the self-energy and for these 
currents. 

The deviations of the anomalous dimension of the 
scalar and pseudoscalar densities from its A^^ = oo value 
are given for the Feynman gauge in Table Ivl 

Thus, in domain- wall fermions for any finite Ns already 
at the one-loop level the anomalous dimensions of the 
continuum and lattice versions of an operator are not 
the same, and this raises some issues about the correct 
procedure with which one must carry out the matching 
of lattice operators to a continuum scheme in this case. 
In fact, the matching formula |5CL.51j 

(g|0'°*|g) " [-lo^oga +R -R J (62) 

is only valid provided the coefficients of the logarithmic 
terms in 

{q\0"^' \q) = (l + 52 ( - JO log aV + R'^') ) • {q\0'^nq) 

(63) 

and 

{q\0'^\q} = (l + .g|^(-7olog^ + i?^)) •(g|0*'--|g) 

(64) 



are the same. But this does not happen for domain-wall 
fermions at finite Ng , because the 70 in Eq. H63|l must be 
replaced with a 7o°'^^ which depends on Ng, and thus 
this logarithmic term cannot be combined with the one 
of Eq. I|64|) , whose coefficient is given by the continuum 
theory. In particular, the p dependence cannot be elimi- 
nated from the matching formula. 

It could well be that all this is connected to the fact 
that the domain-wall theory at finite A^^ is loaded with 
some pathologies, as it possesses no analytical Atiyah- 
Singer index. It is indeed well known that in all formula- 
tions of chiral fermions that recover the correct quantum 
anomalies one needs (in one way or another) an infinite 
number of fermion fields. With a finite number of these 
fields what happens is that either the two chiralities can- 
not be completely separated or that the right anomalies 
cannot be reproduced. In domain- wall at finite Ng, which 
corresponds to a finite number of fermions fields, we could 
then anticipate something like this to happen. We will 
see in the following Section that in the theory truncated 
at finite Ng this is not the only strange feature which 
arises, but that gauge invariance is lost as well, even for 
finite quantities. 



V. RESIDUAL MASS 

In this Section we report, for several choices of Ng (and 
M), the results that we have obtained for Eg. This quan- 
tity enters into the description of the one-loop quark self- 
energy, Eq. H47|l . and determines the critical (or residual) 
mass at this level. The numbers that we have obtained 
for Eg are valid both in the quenched and unquenched 
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TABLE I: Residual mass at tree level in lattice units (multiplied for IGtt^). 



M 


iVs = 8 


Ns = 12 


Ns = 16 


Ns = 20 


Ns = 24 


Ns = 28 


Ns = 32 


Ns = 48 


Ns = oo 


A 1 
U.i 






-0.DD\)(O 


-O.04:( (4: 






1 AOAOO 

-l.UoUz^ 


A 1 AAAA 

-u.iyuyu 


A 



A O 
U.Z 


-\).0Oibi 


-o.yUbbo 


-I.dUUIo 


A K K /I O 


A O A a 


A 1 AAAf! 

-u.iuyyb 


A A/1 KA/1 

-U.U4oU4 


A AA1 OT 

-U.UUlz ( 


A 



0.3 


-4.D4274 


-1.11472 


-U.2d7d4 




A A 1 C /t O 

-U.U1543 


A AAOTA 

-U.UU^yU 


A AAAOA 

-U.UUU89 


A AAAAA 

U.UUUUU 





A A 
0.4 


-i.Dy /oU 


A OOAAA 

-U.z^UUU 


A AOO K 1 


A AAOTA 


A AAA /I O 
-U.UUU40 


A AAAAi? 

-U.UUUUb 


A AAAA1 

-U.UUUUl 


A AAAAA 

U.UUUUU 


A 



0.5 


-U.4dzd4 


A AOOA1 

-0.02891 


A AA 1 O 1 

-0.00181 


A AAA1 1 

-0.00011 


A AAAA1 

-0.00001 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 





0.6 


-U.(J8d9o 


A AAOOO 


A AAAA£? 

-(J.UUUUd 


A AAAAA 

O.UUUUU 


A AAAAA 

U.UUUUU 


A AAAAA 

U.UUUUU 


A AAAAA 

U.UUUUU 


A AAAAA 

U.UUUUU 





A T 
0. ( 


A AAA A O 


A AAAAO 

-U.UUUUo 


A AAAAA 

u.uuuuu 


A AAAAA 

U.UUUUU 


A AAAAA 

U.UUUUU 


A AAAAA 

U.UUUUU 


A AAAAA 

U.UUUUU 


A AAAAA 

U.UUUUU 


A 



0.8 


A AAAOA 

-0.uUu39 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

O.UUUUU 





A A 

o.y 


A AAAAA 
O.UUUUU 


A AAAAA 

O.UUUUU 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A 



1 A 
i.U 


A AAAAA 
O.UUUUU 


A AAAAA 

O.UUUUU 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A 



l.i 


A AAAAA 
O.UUUUU 


A AAAAA 

O.UUUUU 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A 



1.2 


-0.00039 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.3 


-0.00943 


-0.00008 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


-0.08693 


-0.00223 


-0.00006 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


-0.46264 


-0.02891 


-0.00181 


-0.00011 


-0.00001 


0.00000 


0.00000 


0.00000 





1.6 


-1.69750 


-0.22000 


-0.02851 


-0.00370 


-0.00048 


-0.00006 


-0.00001 


0.00000 





1.7 


-4.64274 


-1.11472 


-0.26764 


-0.06426 


-0.01543 


-0.00370 


-0.00089 


0.00000 





1.8 


-9.53767 


-3.90663 


-1.60015 


-0.65542 


-0.26846 


-0.10996 


-0.04504 


-0.00127 





1.9 


-12.91556 


-8.47390 


-5.55973 


-3.64774 


-2.39328 


-1.57023 


-1.03023 


-0.19090 






TABLE 11: CoefRcient of the logarithmic term for Ei, in Feynman gauge. 



M 


Ns=9. 


Ns = 12 


Ns = 16 


iV, = 20 


iV. = 24 


iV, = 28 


Ns = 32 


Ns =48 


= 00 


0.1 





11387 


0.41183 


0.64719 


0.80302 


0.89547 


0.94650 


0.97333 


0.99859 


1 


0.2 





71677 


0.92503 


0.98280 


0.99634 


0.99926 


0.99985 


0.99997 


1.00000 


1 


0.3 





94927 


0.99550 


0.99965 


0.99997 


1.00000 


1.00000 


1.00000 


1.00000 


1 


0.4 





99443 


0.99986 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


0.5 





99964 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


0.6 





99999 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


0.7 


1 


00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


0.8 


1 


00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


0.9 


1 


00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


1.0 


1 


00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


1.1 


1 


00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


1.2 


1 


00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


1.3 


1 


00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


1.4 





99999 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


1.5 





99964 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


1.6 





99443 


0.99986 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


1.7 





94927 


0.99550 


0.99965 


0.99997 


1.00000 


1.00000 


1.00000 


1.00000 


1 


1.8 





71677 


0.92503 


0.98280 


0.99634 


0.99926 


0.99985 


0.99997 


1.00000 


1 


1.9 





11387 


0.41183 


0.64719 


0.80302 


0.89547 


0.94650 


0.97333 


0.99859 


1 



cases, because at one loop internal quark loops can never 
appear in the diagrams that enter in this as well as in 
the other calculations presented in this paper. These di- 
agrams are standard and well known and can be found 
for example in Ref. j4Ql] . 

The behavior of the tadpole diagrams as Ng and M 
change is particularly interesting, and we think that it is 
useful to include here also the values of the tadpole con- 
tributing to El (although in this work we need only the 
results of the tadpole contributing to Eq). Since no pure 
5-dimensional quark propagators appear in the tadpoles, 
for these diagrams the E^t {p) of Eq. H4t)|l is diagonal in the 



fifth dimension and proportional to [i^ ~ Ar/a) G^^{k). 
This is the same integrand of the tadpoles for Wilson 
fermions, and for that action in the case of the tadpole 
diagram contributing to Ei it gave the result (in a general 
covariant gauge) 

Tj = 871^^0 (1-1/4(1 -a)), (65) 

where Zq = 0.154933390231 ... is a well-known integral 
[sif . It is then clear that for domain- wall fermions (where 
we have now r = — 1) the behavior of the tadpole dia- 
grams as a function of Ng (and M) is completely de- 
termined by the damping factors in the fifth dimension. 
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TABLE III: Coefficient of the logarithmic term of the vertex diagram of the vector current, in Feynman gauge. 



A I 


-\',s = 8 


A' 


s = 12 


A' 


s = i() 


A'. = 20 


A" 


, = 21 


A', = 28 


A\ = ;i2 


A' 


. = 18 


A', = -X) 


0.1 


-0 847Q9 


_l 


09908 


_l 


05864 


-1 04859 


_l 


03113 


-1.01779 


-1 00959 




nnn58 

uuuuo 


_\ 


0.2 


-1.04275 


_l 


02123 


_l 


00599 
uuu 


-1 001 H8 

-L . UU -L OO 


_l 


onnsn 


-1 onnn6 

A- . UUUUU 


-1 onnm 

± . UUUU -L 


_x 


nnnnn 

UUUUU 


.\ 


0.3 


-1.01354 


_l 


00158 


_l 


00014 


-1 nnnni 

-L . UUUU -L 


_l 


nnnnn 


-1 nnnnn 

A- . UUUUU 


-1 nnnnn 

J- . UUUUU 


_\ 


nnnnn 

UUUUU 




0.4 


-1.00176 


.1 


00005 


_1 


nnnnn 

UUUUU 


-1 nnnnn 

± .UUUUU 


_x 


nnnnn 

UUUUU 


-1 nnnnn 

± .UUUUU 


-1 nnnnn 

J. .UUUUU 


.\ 


nnnnn 

UUUUU 


_x 


0.5 


-1 nnni9 


.1 


nnnnn 


.1 


nnnnn 


-1 nnnnn 


_1 


nnnnn 


-1 nnnnn 


-1 nnnnn 


.\ 


nnnnn 


.\ 




1 nnnnn 


_1 


00000 

UUUUU 


_1 


00000 
uuuuu 


1 00000 
- -L .uuuuu 


_l 


nnnnn 


1 nnnnn 

- -L .UVJUUU 


1 nnnnn 


_\ 


nnnnn 


_\ 


n 7 


1 nnnnn 


_\ 


00000 
uuuuu 


_\ 


00000 
uuuuu 


1 00000 

- -L .uuuuu 


_l 


nnnnn 

uuuuu 


1 nnnnn 

- -L .uuuuu 


1 nnnnn 

- J- .uuuuu 


_\ 


nnnnn 

uuuuu 


_\ 


u.o 


1 nnnnn 


_\ 


00000 
uuuuu 


_\ 


00000 
uuuuu 


1 00000 

- -L .uuuuu 


_l 


nnnnn 

uuuuu 


1 nnnnn 

- ± .uuuuu 


1 nnnnn 

- ± .uuuuu 


_\ 


nnnnn 

uuuuu 


.\ 


0.9 


-1 nnnnn 


_l 


00000 


_l 


00000 
uuuuu 


-1 00000 

-L .UUUUU 


_l 


nnnnn 

uuuuu 


-1 nnnnn 

± . uuuuu 


-1 nnnnn 

± . uuuuu 


_\ 


nnnnn 

uuuuu 


.\ 


1.0 


-1 nnnnn 


_l 


00000 

uuuuu 


_l 


00000 

uuuuu 


-1 00000 

-L .UUUUU 


_l 


nnnnn 

UUUUU 


-1 nnnnn 

1. . UUUUU 


-1 nnnnn 

L . UUUUU 


_\ 


nnnnn 

UUUUU 


.X 


1.1 


-1 nnnnn 


.1 


nnnnn 

UUUUU 


.1 


nnnnn 

UUUUU 


-1 nnnnn 

J. .UUUUU 


_x 


nnnnn 

UUUUU 


-1 nnnnn 

J. .UUUUU 


-1 nnnnn 

J. .UUUUU 


.\ 


nnnnn 

UUUUU 


.X 




1 nnnnn 

-1. uuuuu 


1 

-1 


nnnnn 
uuuuu 


1 

-1 


nnnnn 
uuuuu 


1 nnnnn 
-1. uuuuu 


1 

-± 


nnnnn 
uuuuu 


1 nnnnn 
-1. uuuuu 


1 nnnnn 
-1. uuuuu 


1 

-1 


nnnnn 
uuuuu 


1 

-X 


1.3 


-1.00000 


-1 


00000 


-1 


00000 


-1.00000 


-1 


00000 


-1.00000 


-1.00000 


-1 


00000 


-1 


1.4 


-1.00000 


-1 


00000 


-1 


00000 


-1.00000 


-1 


00000 


-1.00000 


-1.00000 


-1 


00000 


-1 


1.5 


-1.00012 


-1 


00000 


-1 


00000 


-1.00000 


-1 


00000 


-1.00000 


-1.00000 


-1 


00000 


-1 


1.6 


-1.00176 


-1 


00005 


-1 


00000 


-1.00000 


-1 


00000 


-1.00000 


-1.00000 


-1 


00000 


-1 


1.7 


-1.01354 


-1 


00158 


-1 


00014 


-1.00001 


-1 


00000 


-1.00000 


-1.00000 


-1 


00000 


-1 


1.8 


-1.04275 


-1 


02123 


-1 


00592 


-1.00138 


-1 


00030 


-1.00006 


-1.00001 


-1 


00000 


-1 


1.9 


-0.84792 


-1 


02208 


-1 


05864 


-1.04852 


-1 


03113 


-1.01779 


-1.00952 


-1 


00058 


-1 



TABLE IV: Anomalous dimension of the vector current, in Feynman gauge. 



M 


iV. =8 


iV, =12 


N, = 16 


N, = 20 


AT, = 24 


= 28 


Ns = 32 


N, =48 


A^s = 00 


0.1 


0.73405 


0.61025 


0.41144 


0.24550 


0.13567 


0.07128 


0.03619 


0.00199 





0.2 


0.32598 


0.09620 


0.02312 


0.00504 


0.00104 


0.00021 


0.00004 


0.00000 





0.3 


0.06427 


0.00608 


0.00049 


0.00004 


0.00000 


0.00000 


0.00000 


0.00000 





0.4 


0.00733 


0.00020 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.5 


0.00048 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.6 


0.00002 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.7 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.8 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.9 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.0 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.1 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.2 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.3 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


0.00002 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


0.00048 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.6 


0.00733 


0.00020 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.7 


0.06427 


0.00608 


0.00049 


0.00004 


0.00000 


0.00000 


0.00000 


0.00000 





1.8 


0.32598 


0.09620 


0.02312 


0.00504 


0.00104 


0.00021 


0.00004 


0.00000 





1.9 


0.73405 


0.61025 


0.41144 


0.24550 


0.13567 


0.07128 


0.03619 


0.00199 






Their general effect can be seen by looking at their lead- 
ing contributions for large Ng, which enter the game in 
the combinations 



and 



E 



W, 



Ns-s 



Ns 



2Ns 



s=l 



(66) 



(67) 



These are indeed the leading expressions, in units of 



(68) 



for the tadpole contributions to Eq and Si respectively, 
in the limit of large A''^ . From these expressions one can 
immediately see that the tadpole of Eq vanishes when 
Ns = OO, while the tadpole of Si gives in this limit 
the well-known Wilson number, T;. The damping fac- 
tors thus play a primary role in determining the results 
of the domain- wall tadpoles. We have calculated their 
exact expressions including all subleading terms in Ng, 
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TABLE V: Anomalous dimension of the scalar density, in Feynman gauge. 



M 


iVs =8 


Ns = 12 


Ns = 16 


iV. = 20 


Ns = 24 


iV, = 28 


Ns = 32 


Ns=A8 


Ns= oo 


0.1 


3 27780 


3 67649 


3 58735 


3 391 05 

<J . tji^ i \J' J 


3 22907 


3.12464 


3.06474 


3 00372 


3 


0.2 


3.45424 


3 1 599n 


3 04087 


3 00918 


3 00193 


3 nnn39 


3 nnnns 


3 00000 


3 


0.3 


3 10489 


3 ninsn 


3 00089 


3 00007 


3 00000 


3 nnnnn 


3 nnnnn 


3 00000 


3 


0.4 


3 01260 


3 nnn36 


3 00001 


3 00000 


3 00000 


3 nnnnn 


3 nnnnn 


3 00000 


3 


0.5 


3 00085 


3 onnni 


3 00000 


3 00000 


3 00000 


3 nnnnn 


3 nnnnn 


3 00000 


3 


0.6 


3 00003 


3 nnnnn 

o .KJyjyjyjyj 


3 00000 


3 00000 

o .uuuuu 


3 00000 

O .UUUUU 


3 nnnnn 

O.UUUUU 


3 nnnnn 

O.UUUUU 


3 00000 

O .UUUUU 


3 


0.7 


3 oooon 


3 nnnnn 


3 00000 


3 00000 

O .UUUUU 


3 00000 

O .UUUUU 


3 nnnnn 

O.UUUUU 


3 nnnnn 

O.UUUUU 


3 00000 

O .UUUUU 


3 


0.8 


3 oooon 


3 nnnnn 


3 00000 


3 00000 


3 00000 


3 nnnnn 

O.UUUUU 


3 nnnnn 

O.UUUUU 


3 00000 

O .UUUUU 


3 


0.9 


3 oooon 


3 nnnnn 


3 00000 


3 00000 


3 00000 


3 nnnnn 


3 nnnnn 


3 00000 


3 


1.0 


3 oooon 


3 nnnnn 


3 00000 


3 00000 


3 00000 


3 nnnnn 


3 nnnnn 


3 00000 


3 


1.1 


3 nnnon 


3 nnnnn 


3 00000 


3 00000 


3 00000 


3 nnnnn 


3 nnnnn 


3 00000 


3 


1 9 


o.uuuuu 


Q nnnnn 

o.uuuuu 


Q nnnnn 

o.uuuuu 


Q nnnnn 

o.uuuuu 


Q nnnnn 

o.uuuuu 


Q nnnnn 

o.uuuuu 


Q nnnnn 

o.uuuuu 


Q nnnnn 

O.UUUUU 


■3 

o 


1.3 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3 


1.4 


3.00003 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3 


1.5 


3.00085 


3.00001 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3 


1.6 


3.01260 


3.00036 


3.00001 


3.00000 


3.00000 


3.00000 


3.00000 


3.00000 


3 


1.7 


3.10489 


3.01080 


3.00089 


3.00007 


3.00000 


3.00000 


3.00000 


3.00000 


3 


1.8 


3.45424 


3.15990 


3.04087 


3.00918 


3.00193 


3.00039 


3.00008 


3.00000 


3 


1.9 


3.27780 


3.67649 


3.58735 


3.39105 


3.22907 


3.12464 


3.06474 


3.00372 


3 



TABLE VI: Tadpole contribution to Eq, Eq. 1691 . in Feynman gauge (multiplied for 167r^). 



M 


Ns=8 


Ns = 12 


Ns = 16 


iV, = 20 


Ns = 24 


Ns = 28 


Ns = 32 


Ns =48 


Ns = oo 


0.1 


15.10361 


16.98204 


16.85601 


15.31833 


13.06260 


10.63233 


8.36032 


2.59484 





0.2 


18.33890 


12.98046 


7.73201 


4.17637 


2.12604 


1.04115 


0.49623 


0.02182 





0.3 


12.61592 


4.97585 


1.66803 


0.51427 


0.15080 


0.04277 


0.01184 


0.00006 





0.4 


6.03146 


1.23559 


0.21903 


0.03602 


0.00566 


0.00086 


0.00013 


0.00000 





0.5 


2.10263 


0.20309 


0.01717 


0.00135 


0.00010 


0.00001 


0.00000 


0.00000 





0.6 


0.51309 


0.02003 


0.00069 


0.00002 


0.00000 


0.00000 


0.00000 


0.00000 





0.7 


0.07598 


0.00093 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.8 


0.00476 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.9 


0.00004 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.0 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.1 


-0.00004 


-0.00000 


-0.00000 


-0.00000 


-0.00000 


-0.00000 


-0.00000 


-0.00000 





1.2 


-0.00476 


-0.00001 


-0.00000 


-0.00000 


-0.00000 


-0.00000 


-0.00000 


-0.00000 





1.3 


-0.07598 


-0.00093 


-0.00001 


-0.00000 


-0.00000 


-0.00000 


-0.00000 


-0.00000 





1.4 


-0.51309 


-0.02003 


-0.00069 


-0.00002 


-0.00000 


-0.00000 


-0.00000 


-0.00000 





1.5 


-2.10263 


-0.20309 


-0.01717 


-0.00135 


-0.00010 


-0.00001 


-0.00000 


-0.00000 





1.6 


-6.03146 


-1.23559 


-0.21903 


-0.03602 


-0.00566 


-0.00086 


-0.00013 


-0.00000 





1.7 


-12.61592 


-4.97585 


-1.66803 


-0.51427 


-0.15080 


-0.04277 


-0.01184 


-0.00006 





1.8 


-18.33890 


-12.98046 


-7.73201 


-4.17637 


-2.12604 


-1.04115 


-0.49623 


-0.02182 





1.9 


-15.10361 


-16.98204 


-16.85601 


-15.31833 


-13.06260 


-10.63233 


-8.36032 


-2.59484 






in the rest of the paper, we also show the corresponding 
values for the limiting case of infinite extent in the fifth 
dimension. In the case of Eq we observe that its tad- 
pole contribution presents wide variations with Ns and 
Af , so that in some cases it turns out to be small while 
in other cases it can be substantially large. The tad- 
pole contributing to Si instead has smaller variations. 
This suggests that some care should be used when talk- 
ing about tadpole dominance in relation to domain-wall 
fermions. In fact, the tadpoles contributing to Sq and 
El even decrease toward zero when M ^ or M — > 2, 
as we will see in Sect. IVIIII 



and the tadpole contribution to Eq turns out to be equal 
to 



1 2Ar, 

iV^l + Wo^ ')wq' -2 Wo -j V 



(69) 

while the tadpole contribution to Ei turns out to be equal 
to 



T, 



(1 + w, 



1 Z 



(70) 



Numerical values of these tadpoles for various choices of 
Ns and M are collected in Tables IVll and IVIll where, as 
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TABLE VII: Tadpole contribution to Ei, Eq. I70II . in Feynman gauge (multiplied for 167r^). 



M 


iV. = 8 


Ns = 12 


Ns = 16 


Ns = 20 


A^. = 24 


Ns = 28 


A^, = 32 


Ns =48 


Ns 


= CO 


A 1 
U.i 


6 


88684 


8 


89791 


10 


28165 


11 


14950 


11 


65586 


11 


93540 


12 


08339 


12 


22491 


1 o 

iz 




A O 
U.Z 


10 


71422 


11 


82437 


12 


13713 


12 


21230 


12 


22878 


12 


23220 


12 


23289 


12 


23305 


1 o 

iz 




0.3 


11 


95987 


12 


20785 


12 


23105 


12 


23290 


12 


23304 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


A A 
0.4 


12 


20239 


12 


23224 


12 


23303 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1 o 

iz 


Z66UO 


0.5 


12 


23104 


12 


23304 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


0.6 


12 


23299 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


A T 
0. ( 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1 o 

iz 


OOOAC 

2o6UO 


0.8 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


A A 

o.y 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1 o 

iz 


266UO 


1 A 
i.U 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1 o 

iz 


OOOAK 
ZOOUO 


i.i 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1 o 

iz 




1.2 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1.3 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1.4 


12 


23299 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1.5 


12 


23104 


12 


23304 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1.6 


12 


20239 


12 


23224 


12 


23303 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1.7 


11 


95987 


12 


20785 


12 


23105 


12 


23290 


12 


23304 


12 


23305 


12 


23305 


12 


23305 


12 


23305 


1.8 


10 


71422 


11 


82437 


12 


13713 


12 


21230 


12 


22878 


12 


23220 


12 


23289 


12 


23305 


12 


23305 


1.9 


6 


88684 


8 


89791 


10 


28165 


11 


14950 


11 


65586 


11 


93540 


12 


08339 


12 


22491 


12 


23305 



TABLE VIIL Coefficient of for the contribution of the half-circle diagram to Eq, in Feynman gauge. 



M 


Ns=8 


Ns = 12 


Ns = 16 


Ns = 20 


Ns = 24 


Ns = 28 


Ns = 32 


Ns =48 


A^s = oo 


0.1 


2.75558 


2.69676 


2.39676 


2.00277 


1.60286 


1.24317 


0.94186 


0.26825 





0.2 


3.48374 


2.22615 


1.24197 


0.64382 


0.31901 


0.15331 


0.07208 


0.00306 





0.3 


2.67348 


1.01173 


0.33194 


0.10111 


0.02942 


0.00830 


0.00227 


0.00001 





0.4 


1.52701 


0.31852 


0.05710 


0.00946 


0.00147 


0.00022 


0.00003 


0.00000 





0.5 


0.68927 


0.07235 


0.00641 


0.00051 


0.00004 


0.00000 


0.00000 


0.00000 





0.6 


0.24943 


0.01188 


0.00046 


0.00002 


0.00000 


0.00000 


0.00000 


0.00000 





0.7 


0.07417 


0.00157 


0.00003 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.8 


0.02085 


0.00029 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.9 


0.00770 


0.00012 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.0 


0.00418 


0.00008 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.1 


0.00286 


0.00007 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.2 


0.00304 


0.00007 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.3 


0.01346 


0.00021 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


0.06983 


0.00270 


0.00009 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


0.24400 


0.02258 


0.00186 


0.00013 


0.00001 


0.00000 


0.00000 


0.00000 





1.6 


0.60517 


0.11252 


0.01902 


0.00304 


0.00045 


0.00006 


0.00001 


0.00000 





1.7 


1.12385 


0.36321 


0.10935 


0.03157 


0.00886 


0.00243 


0.00064 


0.00000 





1.8 


1.56704 


0.77955 


0.36471 


0.16644 


0.07499 


0.03349 


0.01486 


0.00053 





1.9 


1.40741 


1.01047 


0.66666 


0.41011 


0.23700 


0.12837 


0.06403 


-0.00472 








Given the strong dependence on M and Ns, tadpole 
improvement (at least in its more common form) seems 
in general not to be an appropriate tool with regard to 
the residual mass, which is here associated with the addi- 
tive mass renormalization arising from the one-loop self- 
energy, Eq. (|48|) . The contribution of the one- loop di- 
agrams to this mass can be inferred from the numbers 
for Eo presented in Tables |3 and IXII (in Feynman and 
Landau gauge respectively), which take into account the 
results for the half-circle diagram given in Tables IVIIII 
and IIXI Since we have performed the calculations in a 
general covariant gauge, we can express all the quantities 



presented in this paper in the form 

A+{l-a)B, (71) 

where A and A + B provide the answer in Feynman and 
Landau gauge respectively, and _B is a number which re- 
mains the same when using fermion formulations rather 
diverse like domain-wall with an infinite extent of the 
fifth dimension, Wilson or overlap. The values of B in 
the case of the contribution of the half-circle diagram to 
Sg at finite are given in Table Hxl If we now compare 
the numbers shown in Tables H ^-nd IXII we can deduce 
that, even after the result for the tadpole is included, the 
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TABLE IX: Coefficient of for the part proportional to 1 ~ a of the contribution of the half-circle diagram to Eq. 



M 


iVs =8 


Ns = 12 


Ns = 16 


iV, = 20 


Ns = 24 


iV. = 28 


Ns = 32 


Ns =48 


Ns = oo 


n 1 
U.l 




Q VVI 01 


O.o4yo4 


O.ODUZ ( 


q AVI 11 

o.U ^ i M 


z.ozu ^ y 


i.yyoo / 


U.OZDDO 


A 

u 


U.z 


4.UU < oi 


z.yoUiz 


1 7AO 1 A 

i. / yziu 




n c;ai ^;a 

u.ouioy 


n O /I Q Q 
U.Z4DOO 


All QA/I 

U.iioU4 


A nncoq 
U.UUOzo 


A 

u 


U.o 


o Q 1 1 op: 
z.oiizo 


1 1 ACi'iK 

i.i4UoO 


U.OO ( ZD 


n 1 OAOfi 
U.izUzD 


n AQ /I Q 
U.Uoo4o 


n A1 AAC 

U.UiUUo 


A nnoon 
U.UUzoU 


A nnnni 
U.UUUUi 


A 
U 






u.zooyu 


U.UoU / o 


n nnQQO 
U.UUooy 


n nni qq 
U.UUloo 


n nnnon 
U.UUUZU 


A nnnnQ 
U.UUUUo 


A Annnn 
U.UUUUU 


A 

u 


U.O 


n /I 'T'T'/i r 
U.4/ /4o 


U.U4D ( 4 


n AAQAQ 

u.uuoyo 


n AAAQO 

U.UUUoz 


n AAAAO 

U.UUUUz 


n AAAAA 

U.UUUUU 


U.UUUUU 


A nnnnn 
U.UUUUU 


A 

u 


U.D 


nil VQV 


U.UU4DO 


n AAA1 
U.UUUiD 


n AAAA1 

U.UUUUi 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


U.UUUUU 


A nnnnn 
U.UUUUU 


A 

u 


U. ( 


U.Ul / oo 


A nnnoo 
U.UUUzz 


n AAAAA 

U.UUUUU 


n AAAAA 

U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


U.UUUUU 


A nnnnn 
U.UUUUU 


A 

u 


U.o 


U.UUiiz 


A nnnnn 
U.UUUUU 


n AAAAA 

U.UUUUU 


n AAAAA 

U.UUUUU 


n AAAAA 

U.UUUUU 


n nnnnn 
U.UUUUU 


A nnnnn 
U.UUUUU 


A nnnnn 
U.UUUUU 


A 

u 


u.y 


U.UUUUl 


A nnnnn 
U.UUUUU 


n AAAAA 

U.UUUUU 


n AAAAA 

U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


A nnnnn 
U.UUUUU 


A nnnnn 
U.UUUUU 


A 

u 


i.U 


u.uuuuu 


A nnnnn 

U.UUUUU 


n nnnnn 

U.UUUUU 


n nnnnn 

U.UUUUU 


n nnnnn 

U.UUUUU 


n nnnnn 

U.UUUUU 


U.UUUUU 


n nnnnn 

U.UUUUU 


n 
u 


i. i 


n nnnni 
-U.UUUUl 


U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


U.UUUUU 


A nnnnn 
U.UUUUU 


A 

u 


1.2 


-U.uulio 


A r\r\r\r\r\ 

0.00000 


0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 





1.3 


-0.02044 


-0.00025 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


-0.13917 


-0.00539 


-0.00018 


-0.00001 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


-0.57387 


-0.05480 


-0.00461 


-0.00036 


-0.00003 


0.00000 


0.00000 


0.00000 





1.6 


-1.65669 


-0.33389 


-0.05873 


-0.00961 


-0.00150 


-0.00023 


-0.00003 


0.00000 





1.7 


-3.49885 


-1.34763 


-0.44675 


-0.13687 


-0.03997 


-0.01130 


-0.00312 


-0.00002 





1.8 


-5.17527 


-3.54181 


-2.07409 


-1.10946 


-0.56133 


-0.27375 


-0.13008 


-0.00568 





1.9 


-4.39070 


-4.73341 


-4.58402 


-4.10098 


-3.46030 


-2.79565 


-2.18639 


-0.67074 






TABLE X: Coefficient of for the complete result of So, in Feynman gauge. 



M 


Ns=8 


Ns = 12 


Ns = 16 


iV, = 20 


Ns = 24 


iV, = 28 


Ns = 32 


Ns =48 


Ns = oo 


0.1 


17.85919 


19.67880 


19.25277 


17.32109 


14.66546 


11.87550 


9.30218 


2.86309 





0.2 


21.82264 


15.20660 


8.97397 


4.82019 


2.44504 


1.19446 


0.56831 


0.02488 





0.3 


15.28940 


5.98758 


1.99997 


0.61538 


0.18022 


0.05106 


0.01412 


0.00007 





0.4 


7.55847 


1.55410 


0.27613 


0.04548 


0.00713 


0.00108 


0.00016 


0.00000 





0.5 


2.79190 


0.27543 


0.02358 


0.00186 


0.00014 


0.00001 


0.00000 


0.00000 





0.6 


0.76252 


0.03191 


0.00115 


0.00004 


0.00000 


0.00000 


0.00000 


0.00000 





0.7 


0.15015 


0.00250 


0.00004 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.8 


0.02561 


0.00031 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.9 


0.00774 


0.00012 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.0 


0.00418 


0.00008 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.1 


0.00283 


0.00007 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.2 


-0.00172 


0.00006 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.3 


-0.06252 


-0.00073 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


-0.44327 


-0.01733 


-0.00060 


-0.00002 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


-1.85863 


-0.18051 


-0.01531 


-0.00122 


-0.00009 


-0.00001 


0.00000 


0.00000 





1.6 


-5.42629 


-1.12306 


-0.20001 


-0.03298 


-0.00521 


-0.00080 


-0.00012 


0.00000 





1.7 


-11.49207 


-4.61264 


-1.55867 


-0.48269 


-0.14194 


-0.04034 


-0.01120 


-0.00006 





1.8 


-16.77186 


-12.20090 


-7.36730 


-4.00994 


-2.05105 


-1.00765 


-0.48138 


-0.02129 





1.9 


-13.69620 


-15.97157 


-16.18935 


-14.90821 


-12.82561 


-10.50396 


-8.29629 


-2.59956 






contribution to Eg proportional to 1 — a is in general not 
equal to zero, which is what one would have instead ob- 
tained if A^s = oo or if Wilson or overlap fermions were 
used. All this means that at finite Ng the residual mass, 
which is derived from Sq according to Eq. H48|l . is not 
a gauge invariant quantity anymore. Although numeri- 
cally the deviations from gauge invariance remain in most 
cases rather small (because there are large cancellations 
between the contributions of the half-circle and tadpole 
diagrams), we encounter here another of the pathologi- 
cal features of the domain- wall theory truncated at finite 
Ns- Thus, anomalous dimensions as well as terms pro- 



portional to 1 — a, two of the quantities that remain the 
same when using a wide variety of fermionic actions, as- 
sume instead new values when the theory of domain-wall 
fermions is truncated at a finite Ns ■ 

It is also interesting to compare the one-loop results 
presented here for Sq to the numbers that one obtains us- 
ing Wilson fermions, which are —51.43471 for the unim- 
proved and —31.98644 for the improved (with Csw = 1) 
case. We can then see that the domain- wall results are 
much smaller than the Wilson numbers, even for Ns as 
small as 8. 

The values of awrel which come out from our results 
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TABLE XI: Coefficient of for the complete result of Eo, in Landau gauge. 



M 


iVs =8 


Ns = 12 


Ns = 16 


iV, = 20 


Ns = 24 


Ns = 28 


Ns = 32 


Ns =48 


Ns = oo 


n 1 
U.l 


1 / .zoU / U 


ly.zuozu 


1 Q QQQQl 

lo.f5oooi 


1 7 np;i 7fi 
i / .Uol / o 


1/1 /1 71 r^Q 
14.4 ( iOo 


1 1 7QQOn 

ii. / oozU 


o onc;o7 

y.zuoy ( 


Z.o4iUD 


A 

u 


U.z 


zi.z4o / z 


1 A m 1 (^1 
i4.yiiDi 


o.oooUo 


A 7P: /I Q A 

4. /04o4 


O /I 1 c^oo 
Z.4iOZZ 


1 1 Q 1 AA 

i.ioiUU 


A c;^;ooA 
U.ODZZy 


U.Uz4dd 


U 


U.o 


i4.y4DDo 


O.oooy / 


i.y / Uzz 


U.dU rUo 


A 1 77AC^ 

[J.Li / yo 


A ACA/I C 

U.UoU4o 


A A1 QCid 

u.uioyu 


A AAAA7 
U.UUUU / 


U 


n A 


T ACiOQ 1 

/ .4uyoi 


i.oZyiU 


U.Z ( Zlo 


n O/l /I Q7 
U.U44o / 


n nn7n/i 
U.UU /U4 


n nni n7 
U.UUlU / 


n nnni « 
U.UUUiD 


A AAAAA 

U.UUUUU 


A 

U 


U.O 


z. ^ 4oDy 


n 071 /1 1 
U.z { i4i 


n A0Q07 

U.Uzoz / 


A AAI QA 

U.UUio4 


A AAA1 A 

U.UUUi4 


A AAAA1 
U.UUUUi 


A AAAAA 
U.UUUUU 


A AAAAA 
U.UUUUU 


U 


U.D 


n 7c; 1 i^o 


U.Uoloo 


U.UUi 14 


U.UUUU4 


n nnnnn 
U.UUUUU 


n f\c\f\c\f\ 
U.UUUUU 


n nnnnn 
U.UUUUU 


n c\f\c\f\c\ 
U.UUUUU 


A 

U 


U. i 


n 1 /I Q71 


u.uuz4y 


U.UUUU4 


n c\f\c\f\c\ 
U.UUUUU 


n nnnnn 
U.UUUUU 


n f\c\f\c\f\ 
U.UUUUU 


n nnnnn 
U.UUUUU 


n c\f\c\f\c\ 
U.UUUUU 


A 

U 


U.o 


U.Uzoo4 


U.UUUol 


U.UUUU i 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


A AAAAA 

U.UUUUU 


A 

u 


u.y 


U.UU M 4 


U.UUUIZ 


U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


n nnnnn 
U.UUUUU 


A AAAAA 

U.UUUUU 


A 

u 


i.U 


U.UU410 


U.UUUUo 


n nnnnn 

U.UUUUU 


A r\r\r\r\r\ 
U.UUUUU 


n nnnnn 

U.UUUUU 


A r\r\r\r\r\ 
U.UUUUU 


A nnnnn 

U.UUUUU 


A r\r\r\r\r\ 
U.UUUUU 


A 

u 


i. i 


U.UUZoo 


U.UUUU / 


n nnnnn 
U.UUUUU 


n C\f\c\f\c\ 
U.UUUUU 


n nnnnn 
U.UUUUU 


n f\c\f\c\f\ 
U.UUUUU 


n nnnnn 
U.UUUUU 


n C\f\c\f\c\ 
U.UUUUU 


A 

u 


1.2 


-U.U0179 


O.OOOOd 


0.00000 


0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 


A AAAAA 

0.00000 





1.3 


-0.06397 


-0.00074 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


-0.45417 


-0.01771 


-0.00061 


-0.00002 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


-1.90684 


-0.18453 


-0.01562 


-0.00124 


-0.00009 


-0.00001 


0.00000 


0.00000 





1.6 


-5.57512 


-1.14806 


-0.20399 


-0.03359 


-0.00530 


-0.00081 


-0.00012 


0.00000 





1.7 


-11.83694 


-4.71631 


-1.58842 


-0.49100 


-0.14421 


-0.04095 


-0.01136 


-0.00006 





1.8 


-17.36240 


-12.49760 


-7.50839 


-4.07530 


-2.08087 


-1.02111 


-0.48740 


-0.02151 





1.9 


-14.31100 


-16.45947 


-16.55937 


-15.17962 


-13.02025 


-10.64153 


-8.39260 


-2.62158 






TABLE XII: Residual mass in lattice units at /3 = 6 in Feynman gauge. 



M 


Ns=8 


Ns = 12 


Ns = 16 


iV, = 20 


Ns = 24 


iV, = 28 


Ns = 32 


iV. =48 


Ns = oo 


0.1 


-0.23258 


-0.21982 


-0.19777 


-0.16935 


-0.13898 


-0.11021 


-0.08507 


-0.02538 





0.2 


-0.24466 


-0.15313 


-0.08590 


-0.04485 


-0.02234 


-0.01078 


-0.00508 


-0.00022 





0.3 


-0.15850 


-0.05761 


-0.01858 


-0.00560 


-0.00162 


-0.00045 


-0.00012 


0.00000 





0.4 


-0.07457 


-0.01452 


-0.00251 


-0.00041 


-0.00006 


-0.00001 


0.00000 


0.00000 





0.5 


-0.02650 


-0.00251 


-0.00021 


-0.00002 


0.00000 


0.00000 


0.00000 


0.00000 





0.6 


-0.00699 


-0.00028 


-0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.7 


-0.00133 


-0.00002 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.8 


-0.00022 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.9 


-0.00007 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.0 


-0.00004 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.1 


-0.00002 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.2 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.3 


0.00047 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


0.00319 


0.00013 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


0.01276 


0.00134 


0.00012 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 





1.6 


0.03507 


0.00809 


0.00151 


0.00026 


0.00004 


0.00001 


0.00000 


0.00000 





1.7 


0.06763 


0.03189 


0.01147 


0.00367 


0.00110 


0.00032 


0.00009 


0.00000 





1.8 


0.08121 


0.07828 


0.05207 


0.02971 


0.01562 


0.00781 


0.00378 


0.00017 





1.9 


0.03385 


0.08119 


0.10149 


0.10278 


0.09314 


0.07875 


0.06353 


0.02074 






for Eo are reported in Tables IXIll to IXVl One can imme- 
diately notice that the deviations from the case of exact 
chiral symmetry are much more pronounced when M is 
close to 0.1 or 1.9. Since M = 1.8 is at present a pre- 
ferred choice for Monte Carlo simulations, as it appears 
to minimize chiral violations at the nonperturbative level 
01 , we can focus on this value of M and observe that 
at one loop the critical mass seems to be still large up to 
Ng = 16 or Ns = 20, becoming then somewhat smaller 
for Ns — 24 and higher. Since the tree level and one-loop 
contributions are of a different sign, it could also be that 
the two-loop expression enters again with a negative sign 



and that its effect is to damp the chiral violations that 
we have obtained. But if higher-order corrections do not 
at the end strongly compensate the results of the resid- 
ual mass at one loop, all this suggests that values like 
Ns = 24 would be better choices for the simulations. 

The region around M = 1.8 suffers perturbatively 
rather badly from another problem as well. When M is 
indeed close to 0.1 or 1.9, the numerical convergence of 
the integrals (which are computationally quite demand- 
ing, being in this case sums over a six-dimensional space) 
is much worse then when M is nearer to 1. This behav- 
ior had already been noticed in the calculations for the 



16 



TABLE XIII: Residual mass in lattice units at /3 = 6 in Landau gauge. 



M 


iVs =8 


Ns = 12 


Ns = 16 


iV, = 20 


Ns = 24 


iV. = 28 


Ns = 32 


iV. =48 


Ns = oo 


0.1 


-0 22770 


-0 21 582 


-0 1 9469 


-0 16708 


-0 1 3735 


-0 1 0905 


-0 08425 


-0 02520 

1^. yj^fj^yj 





0.2 


-0 23978 


-0 15064 


-0.08471 


-0 04430 


-0 02209 


-0 01067 


-0 00503 

yj • yjyjtjyj'j 


-0 00022 

1^. yjyjyj^^ 





0.3 


-0 1 5560 


-0 05674 


-0 01 833 


-0 00553 


-0 00160 

Kj • yjyj -L vjvv 


-0 00045 

yj • yjyjyj^tj 


-0 00012 

yj • yjyjyj -Lij 


00000 

yj • yjyjyjyjyj 





0.4 


-0 07331 


-0 01430 


-0 00248 

\J • \J\J iJ^Vj 


-0 00040 


-0 00006 

yj • yjyjyjyjyj 


-0 00001 

yj • yjyjyjyj ±. 


00000 

yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 





0.5 


-0 02610 


-0 00247 

\J • \J\J 1 


-0 00021 


-0 00002 


00000 


00000 

yj . yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 





0.6 


-0 nnfi9n 


-0 00028 


-0 00001 


00000 


00000 

yj •\J\J\JU\J 


00000 
VJ .uuuuu 


00000 

U.UUUUU 


00000 

U.UUUUU 





0.7 


-0 nm "V) 


-0 00002 


00000 


00000 

yj • UUKJKJXJ 


00000 

yj • UKJKJKJKJ 


00000 
VJ . uuuuu 


00000 

U.UUUUU 


00000 
u . yjyjyjyjyj 





0.8 


-0 00022 


00000 


00000 


00000 

yj • UKJKJKJXJ 


00000 

yj • U\J\J\J\J 


00000 

yj • uuuuu 


00000 

U.UUUUU 


00000 
yj .yjyjyjyjyj 





0.9 


-0 00007 


00000 


00000 


00000 


00000 

yj . L/ L/ 


00000 

yj . yjyjyj'^yj 


00000 

yj • yjyjyjyjyj 


00000 

yj . yjyjyjyjyj 





1.0 


-0 00004 


00000 


00000 


00000 


00000 


00000 

yj . yjyjyj'^yj 


00000 

yj • yjyjyjyjyj 


00000 

yj . yjyjyjyjyj 





1.1 


-0 00002 


00000 


00000 


00000 


00000 


00000 

yj . yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 





1 9 


n nnnni 

U.UUUUl 




n nnnnn 


n nnnnn 


n nnnnn 

U.UUUUU 


n nnnnn 

U.UUUUU 


n nnnnn 

U.UUUUU 


n nnnnn 

U.UUUUU 


n 
u 


1.3 


0.00048 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


0.00328 


0.00014 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


0.01317 


0.00137 


0.00012 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 





1.6 


0.03632 


0.00830 


0.00154 


0.00026 


0.00004 


0.00001 


0.00000 


0.00000 





1.7 


0.07054 


0.03276 


0.01172 


0.00374 


0.00112 


0.00032 


0.00009 


0.00000 





1.8 


0.08620 


0.08078 


0.05326 


0.03026 


0.01587 


0.00793 


0.00383 


0.00017 





1.9 


0.03905 


0.08531 


0.10461 


0.10507 


0.09478 


0.07991 


0.06434 


0.02093 






TABLE XIV: Residual mass in lattice units at /3 = 5.2 in Feynman gauge. 



M 


Ns=8 


Ns = 12 


Ns = 16 


Ns = 20 


Ns = 24 


Ns = 28 


Ns = 32 


Af. =48 


As = oo 


0.1 


-0.25578 


-0.24538 


-0.22278 


-0.19185 


-0.15803 


-0.12564 


-0.09715 


-0.02910 





0.2 


-0.27300 


-0.17289 


-0.09756 


-0.05111 


-0.02552 


-0.01233 


-0.00582 


-0.00025 





0.3 


-0.17836 


-0.06539 


-0.02118 


-0.00640 


-0.00185 


-0.00052 


-0.00014 


0.00000 





0.4 


-0.08439 


-0.01653 


-0.00287 


-0.00047 


-0.00007 


-0.00001 


0.00000 


0.00000 





0.5 


-0.03013 


-0.00287 


-0.00024 


-0.00002 


0.00000 


0.00000 


0.00000 


0.00000 





0.6 


-0.00798 


-0.00033 


-0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.7 


-0.00152 


-0.00002 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.8 


-0.00025 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





0.9 


-0.00008 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.0 


-0.00004 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.1 


-0.00003 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.2 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.3 


0.00055 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


0.00377 


0.00015 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


0.01518 


0.00158 


0.00014 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 





1.6 


0.04212 


0.00955 


0.00177 


0.00030 


0.00005 


0.00001 


0.00000 


0.00000 





1.7 


0.08256 


0.03788 


0.01349 


0.00430 


0.00129 


0.00037 


0.00010 


0.00000 





1.8 


0.10300 


0.09413 


0.06164 


0.03492 


0.01828 


0.00912 


0.00440 


0.00020 





1.9 


0.05165 


0.10194 


0.12252 


0.12214 


0.10980 


0.09239 


0.07430 


0.02412 






Ns = oo limit presented in j4flj , however now the conver- 
gence is even slower than in that case, and it also tends 
to become worse when Ns is small. 

The values of the residual mass amrel , that is after the 
one-loop corrections are added to the tree-level expres- 
sion, have the same sign of the critical mass of Wilson 
fermions only for M > 1.2. In this region of M they are 
positive, at least for even Ns and if the coupling is not 
very small. If one looks at the columns corresponding to 
Ns = 8, our results seem to indicate that the minimal 
amount of chiral violations is attained for M ~ 1.2 (in 
Feynman or Landau gauge). The optimal choice of M 



from the point of view of 1-loop calculations would then 
be around these values. This effect can be related to the 
renormalization of M, a quantity which is not protected 
by chiral symmetry and is moved by radiative corrections 
away from its free value, M — 1 (where the tree-level 
residual mass vanishes). One can conjecture that higher- 
loop corrections and nonperturbative effects would shift 
this optimal value further on, until the chiral violations 
approach a minimal point around M — 1.8, which seems 
to be the choice that provides the smallest residual mass 
in Monte Carlo simulations. It could be that the results 
of the one-loop calculations for M = 1.8 are not substan- 
tially changed once higher-loop corrections are included. 
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TABLE XV: Residual mass in lattice units at /3 = 5.2 in Landau gauge. 



M 


iVs =8 


Ns = 12 


Ns = 16 


iV, = 20 


Ns = 24 


iV. = 28 


Ns = 32 


iV. =48 


Ns = oo 


0.1 


-0 25014 


-0.24077 


-0 21 923 

yj m ^ A- fj 


-0 18923 


-0.15614 


-0 1 2430 


-0 09621 

yj • yj ijyjij ±. 


-0 02889 





0.2 


-0 26738 


-0.17001 


-0 0961 9 

yj .yj rjyj 1- rj 


-0 05047 

yj .yj'jyj^ ■ 


-0 02523 

yj . yj£jrj^tj 


-0 01220 

yj ^yj A-^^yj 


-0 00576 

yj • yjyj\j i v.^ 


-0 00025 

yj . yjyjyj £j\j 





0.3 


-0 1 7502 


-0 06438 


-0 02089 

yj . yj^yjUiJ 


-0 00632 

\J .yjyjyjtj^ 


-0 00183 

yj • yjyj i-yjtj 


-0 00051 

yj • yjyjyj\j ±. 


-0 00014 

yj • yjyjyj -Lt: 


00000 

yj . yjyjyjyjyj 





0.4 


-0 08294 


-0 01 629 


-0 00283 


-0 00046 

yj . yjyjyj^yj 


-0 00007 

yj • yjyjyjyj i 


-0 00001 

yj • yjyjyjyj ±. 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 





0.5 


-0 02966 


-0 00283 


-0 00024 

v7 . yjyjyj^^ 


-0 00002 

v7 . yjyjyjyj^ 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj . yjyjyjyjyj 





0.6 


-0 00787 
\j .yjyj 1 (J 1 


-0 000S2 


-0 00001 


00000 


00000 

yJ •yjyjyjyjyj 


00000 

yj •yjyjyjyjyj 


00000 

U. VJUVJUU 


00000 
u.uuuuu 





0.7 


-0 001 51 


-0 00002 


00000 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj .yjyjyjyjyj 


n 00000 

U . yjyjyjyjyj 





0.8 


-0 00025 


00000 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj .yjyjyjyjyj 


00000 

yj •yjyjyjyjyj 





0.9 


-0 00008 


00000 


00000 

yj . yjyjyjyjyj 


00000 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj . yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 





1.0 


-0 00004 


00000 


00000 

. yjyjyjyjyj 


00000 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj . yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 





1.1 


-0 00003 


00000 


00000 

. yjyjyjyjyj 


00000 


00000 

yj • yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 


00000 

yj . yjyjyjyjyj 


00000 

yj • yjyjyjyjyj 





1 9 


n 00009 


00000 


00000 


00000 


00000 


00000 
u.uuuuu 


00000 

U.UUUUU 


00000 
u.uuuuu 


n 
u 


1.3 


0.00056 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.4 


0.00387 


0.00016 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 





1.5 


0.01565 


0.00161 


0.00014 


0.00001 


0.00000 


0.00000 


0.00000 


0.00000 





1.6 


0.04357 


0.00979 


0.00181 


0.00030 


0.00005 


0.00001 


0.00000 


0.00000 





1.7 


0.08592 


0.03889 


0.01378 


0.00438 


0.00131 


0.00038 


0.00011 


0.00000 





1.8 


0.10875 


0.09702 


0.06302 


0.03555 


0.01857 


0.00925 


0.00446 


0.00020 





1.9 


0.05763 


0.10669 


0.12612 


0.12479 


0.11169 


0.09373 


0.07524 


0.02433 






TABLE XVI: Residual mass at (3 ■ 
Landau gauge. 



6.0 for some odd Ns, in 



TABLE XVII: Residual mass at /? : 
Landau gauge. 



5.2 for some odd Ns, in 



M 


iV. = 13 


iV. = 14 iV, = 15 


Ns = 16 Af, = 17 Ns 


= 18 


0.1 


-0.211 


-0.206 


-0.201 


-0.195 


-0.188 


-0.181 


0.2 


-0.132 


-0.114 


-0.099 


-0.085 


-0.072 


-0.062 


0.3 


-0.043 


-0.033 


-0.025 


-0.018 


-0.014 


-0.010 


0.4 


-0.009 


-0.006 


-0.004 


-0.002 


-0.002 


-0.001 


0.5 


-0.001 


-0.001 


0.000 


0.000 


0.000 


0.000 


0.6 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.7 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.8 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.9 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.0 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.1 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.2 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.3 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.4 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.5 


-0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


1.6 


-0.006 


0.004 


-0.002 


0.002 


-0.001 


0.001 


1.7 


-0.026 


0.020 


-0.015 


0.012 


-0.009 


0.007 


1.8 


-0.074 


0.067 


-0.060 


0.053 


-0.047 


0.041 


1.9 


-0.092 


0.098 


-0.102 


0.105 


-0.106 


0.107 



M 


Ns = 13 


Ns = 14 Af, = 15 


Ns ^IQ Ns = 17 Ns 


= 18 


0.1 


-0.236 


-0.231 


-0.226 


-0.219 


-0.212 


-0.205 


0.2 


-0.149 


-0.129 


-0.112 


-0.096 


-0.082 


-0.070 


0.3 


-0.049 


-0.037 


-0.028 


-0.021 


-0.016 


-0.012 


0.4 


-0.011 


-0.007 


-0.004 


-0.003 


-0.002 


-0.001 


0.5 


-0.002 


-0.001 


0.000 


0.000 


0.000 


0.000 


0.6 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.7 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.8 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.9 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.0 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.1 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.2 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.3 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.4 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


1.5 


-0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


1.6 


-0.006 


0.004 


-0.003 


0.002 


-0.001 


0.001 


1.7 


-0.030 


0.024 


-0.018 


0.014 


-0.010 


0.008 


1.8 


-0.089 


0.080 


-0.071 


0.063 


-0.055 


0.048 


1.9 


-0.114 


0.120 


-0.124 


0.126 


-0.127 


0.127 



while for M < 1.8 they are increased by renormahza- 
tion so that the minimum finally ends up located near 
M = 1.8. 

Notice that for M = 1.9 the residual mass at Ng — 12 
is larger than at Ng — 8, and at Ns = 16 is even larger. 
We will discuss more in detail this kind of behavior in 

Sect, rvml 

As an aside, we have also investigated what happens 
when one chooses an odd value of Ng- As can be seen 
from Tables IXVII and IXVIII which are representative of 
the general situation, when M is smaller than 1 the resid- 



ual mass has always a negative sign for Ng even or odd, 
but in the more interesting case of M greater than 1 the 
one-loop residual mass turns out to be positive when Ng 
is even and negative when Ng is odd. This could be re- 
lated to the su gges tion which was made by Neuberger 
some years ago |53|, according to which for odd Ng the 
theory which is simulated could correspond to the 6 = tt 
regime of QCD. 

The residual mass obviously changes also when the 
coupling go is varied, as an effect of the loop correc- 
tions. Indeed, the overlap between the chiral modes liv- 
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TABLE XVIII: Coefficient of for the quantity A = Zv — Za ~ —{Zs — Zp)/2, which vanishes when chiral symmetry is 
restored at A^'s = co, in Feynman gauge. 



M 


iVs =8 


Ns = 12 


Ns = 16 


Ns = 20 


Ns = 24 


Ns = 28 


Ns = 32 


Ns = 48 


Ns = oo 


n 1 
U.i 


6 


9399 


A A AO A 

4.44o4 


O CiC\A 1 


i.O i 4y 


n QA1 a 
U.oylD 


n A QAA 

u.45yy 


A Ol^OO 

U.zDzz 


A A1 7A 
U.Ui (\} 


U 


U.z 


2 


3272 


n vol 1 


U.iy40 


A A/1 

U.U4j' i 


A A1 1 A 
U.UllU 


A AAOC; 

U.UUzO 


A AAAC; 

U.UUUO 


A AAAA 
U.UUUU 


U 


U.o 





5894 


U.UDDo 


U.UUD4 


A AAAi^ 
U.UUUD 


n AAAA 

U.UUUU 


A AAAA 
U.UUUU 


A AAAA 

U.UUUU 


A AAAA 
U.UUUU 


U 


A A 

U.4 





0994 


U.UUo4 


A nnni 
U.UUUi 


A Annn 
U.UUUU 


n AAAA 

U.UUUU 


A AAAA 
U.UUUU 


A AAAA 

U.UUUU 


A AAAA 
U.UUUU 


U 


U.o 





0106 


U.UUUl 


A nnnn 
U.UUUU 


A nnnn 
U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


A 

U 


U.O 





0007 


U.UUUU 


A AAnn 
U.UUUU 


A Annn 
U.UUUU 


n AAAA 

U.UUUU 


A AAAA 
U.UUUU 


A AAAA 
U.UUUU 


A AAAA 
U.UUUU 


U 


A 7 





0000 


n nnnn 
U.UUUU 


A nnnn 
U.UUUU 


A nnnn 
U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


A 

U 


A Q 
U.O 





0000 


n nnnn 
U.UUUU 


A nnnn 
U.UUUU 


A nnnn 
U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


A 

U 


A O 

u.y 





0000 


n nnnn 
U.UUUU 


A nnnn 
U.UUUU 


A nnnn 
U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


A 

u 


1 A 

i.U 





0000 


n nnnn 
U.UUUU 


A nnnn 
U.UUUU 


A nnnn 
U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


n AAAA 

U.UUUU 


A 

u 


1.1 





0000 


A f\f\f\r\ 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

o.uooo 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 





1.2 





0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.3 





0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.4 





0005 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.5 





0083 


0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.6 





0701 


0.0029 


0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.7 





3140 


0.0501 


0.0056 


0.0005 


0.0000 


0.0000 


0.0000 


0.0000 





1.8 





8354 


0.3672 


0.1288 


0.0374 


0.0096 


0.0022 


0.0005 


0.0000 





1.9 


2 


2075 


1.5323 


1.0348 


0.6904 


0.4490 


0.2811 


0.1686 


0.0153 






TABLE XIX: Coefficient of for the part proportional to 1 — a of A. 



M 


Ns=8 


iV. = 12 


A^. = 16 


Ns = 20 


Af. =24 


Ns = 28 


Ns = 32 


Af, =48 


A^s = oo 


0.1 


-1.1115 


-0.8058 


-0.5284 


-0.3257 


-0.1914 


-0.1081 


-0.0590 


-0.0042 





0.2 


-0.4075 


-0.1438 


-0.0416 


-0.0106 


-0.0025 


-0.0006 


-0.0001 


0.0000 





0.3 


-0.1092 


-0.0140 


-0.0014 


-0.0001 


0.0000 


0.0000 


0.0000 


0.0000 





0.4 


-0.0191 


-0.0007 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.5 


-0.0020 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.6 


-0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.7 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.8 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.9 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.1 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.2 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.3 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.4 


-0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.5 


-0.0022 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.6 


-0.0186 


-0.0008 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.7 


-0.0879 


-0.0134 


-0.0015 


-0.0001 


0.0000 


0.0000 


0.0000 


0.0000 





1.8 


-0.2750 


-0.1064 


-0.0353 


-0.0100 


-0.0025 


-0.0006 


-0.0001 


0.0000 





1.9 


-0.9106 


-0.5678 


-0.3473 


-0.2140 


-0.1314 


-0.0790 


-0.0461 


-0.0040 






ing near the two walls depends also on the strength of 
the gauge coupling, and for strong couplings this overlap 
tends to acquire some nonnegligible values. Furthermore, 
the residual mass explicitly depends also on the value of 
the lattice spacing a, as can be seen from Eq. H48|l . since 
all the numbers presented in the Tables are given, as else- 
where in the paper, in lattice units. The residual mass is 
thus different for quenched and unquenched simulations 
made at the same lattice spacing, because a and go are 
related in a different way. Let us now consider some typi- 
cal values of the parameters at which simulations are cur- 



rently performed. We first fix M = 1.8 and Ng ~ 16. In 
the quenched case we can then take (3 — 6.0, which corre- 
sponds to a lattice spacing of about 2 GeV, and according 
to our one-loop calculations this gives arrires = 0.05207 in 
Feynman gauge and arrires = 0.05326 in Landau gauge, 
that is TUres ~ 104 MeV and m,res ^ 107 MeV respec- 
tively. The dependence on the gauge seems thus from 
the numerical point of view not to be very significant. 

Notice also that if mrel in Eq. ISHIwere positive, the val- 
ues of Tflres would be even larger. In the unquenched 
case, if we take /3 = 5.2, which corresponds roughly to 
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TABLE XX: Coefficient of for the coefficient of the power-divergent mixing of O^j , Cmix in Eq. 17511 . in Feynman gauge. 



M 


iVs = 8 


Ns = 12 


Ns = 16 


iV, = 20 


Ns = 24 


Ns = 28 


Ns = 32 


Ns = 48 


Ns = oo 


A 1 
U.i 


1/1 1 O /I /I 

-i4.iz44 


-o.444o 


-O.ODOZ 


Q AQO A 

-0.4804 


o oo A a 
-z.Zo4d 


1 A f\f\Ci 


A AO A O 


A 1 OOA 


A 



A O 
U.Z 


KOI 


-2.4Ui4 


A AA1 A 

-0.9910 


-0.40/0 


A 1 d'y A 


-U.Uboo 


A AOO Q 

-U.Uzo^ 


A AAAO 

-U.UUUo 


A 



0.3 


O /I TO O 

-2.4733 


-0.6070 


-0.1475 


-0.0359 


-U.UU87 


-0.0021 


-0.0005 


A AAAA 

U.UUUU 





A A 
0.4 


A A O 1 


A 1 O K 

-O.izDO 


A A1 f?A 

-0.01d9 


A AAOO 

-0.0022 


-U.UUUo 


A AAAA 
U.UUUU 


A AAAA 
U.UUUU 


A AAAA 

U.UUUU 


A 



0.5 


-0.3094 


-0.0213 


-0.0014 


-0.0001 


r\ AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 





0.6 


-0.0893 


-0.0029 


-0.0001 


n r\r\r\r\ 

u.uuuu 


A AAAA 

(J.UUUU 


A AAAA 

U.UUUU 


A AAAA 

U.UUUU 


A AAAA 

U.UUUU 





A 

0. ( 


A AOQ 1 

-O.Oz^i 


A AAA /I 

-0.0004 


A AAAA 
0.0000 


U.UUUU 


A AAAA 

U.UUUU 


A AAAA 

U.UUUU 


A AAAA 

U.UUUU 


A AAAA 

U.UUUU 


A 



0.8 


-0.0059 


-0.0001 


A AAAA 

0.0000 


A r\r\r\r\ 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 





A A 

o.y 


A AA1 K 

-O.OOio 


A AAAA 
0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A 



1 A 
1.0 


A AAAA 
0.0000 


A AAAA 
0.0000 


A AAAA 
0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A 



l.i 


A AAAO 

O.OOOo 


A AAAA 
0.0000 


A AAAA 
0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A AAAA 

0.0000 


A 



1.2 


0.0014 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.3 


0.0031 


0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.4 


0.0127 


0.0003 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.5 


0.0546 


0.0029 


0.0002 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.6 


0.1969 


0.0207 


0.0024 


0.0003 


0.0000 


0.0000 


0.0000 


0.0000 





1.7 


0.5774 


0.1129 


0.0229 


0.0049 


0.0011 


0.0002 


0.0001 


0.0000 





1.8 


1.3844 


0.4817 


0.1715 


0.0615 


0.0224 


0.0083 


0.0031 


0.0001 





1.9 


2.9825 


1.5917 


0.9340 


0.5665 


0.3480 


0.2148 


0.1328 


0.0198 






TABLE XXI: Coefficient of for the part proportional to 1 — a of c, 



M 


iV. =8 


iV. = 12 


Ns = 16 


Ns = 20 


Ns=2A 


Ns = 28 


Ns = 32 


iVs =48 


Ns = oo 


0.1 


0.3241 


0.1056 


0.0395 


0.0162 


0.0071 


0.0033 


0.0015 


0.0001 





0.2 


0.0480 


0.0093 


0.0019 


0.0004 


0.0001 


0.0000 


0.0000 


0.0000 





0.3 


0.0093 


0.0008 


0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.4 


0.0014 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.5 


0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.6 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.7 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.8 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





0.9 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.1 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.2 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.3 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.4 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.5 


-0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.6 


-0.0013 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.7 


-0.0060 


-0.0007 


-0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 





1.8 


0.0035 


-0.0043 


-0.0015 


-0.0004 


-0.0001 


0.0000 


0.0000 


0.0000 





1.9 


0.4015 


0.0982 


0.0218 


0.0029 


-0.0011 


-0.0013 


-0.0009 


-0.0001 






the same lattice spacing of 2 GeV but now to a different 
bare coupling, we then obtain arrires — 0.06164 in Feyn- 
man gauge and arrires — 0.06302 in Landau gauge, that is 
rrires ~ 123 MeV and rrires 126 MeV respectively. Re- 
sults from dynamical domain- wall simulations show that 
the residual mass in full QCD assumes much higher val- 
ues than in the quenched case. Dynamical domain- wall 
fermions present then a larger explicit chiral symmetry 
breaking. We can infer a similar effect which goes in the 
same direction from our one-loop perturbative results, 
although here the difference between quenched and full 
QCD is not as pronounced as it has been observed in 



the simulations, suggesting that higher-order corrections 
and especially nonperturbative effects presumably play a 
significant role. 

We have automated the calculations presented in this 
article by developing suitable FORM codes 'Hs'l, which 
are a sizeable extension of the programs which were used 
in the calculations presented in Ref. for Ns = oo. 
These codes are now able to compute matrix elements 
for general values of Ng and M. 

The integrals of the Feynman diagrams have been nu- 
merically evaluated by keeping Ns fixed while at the 
same time refining the integration grid in the usual four- 
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dimensional momentum space (where we needed to con- 
sider values of = 120). One might wonder whether 
this way of performing the computation of the integrals, 
which amounts to taking a finite extra dimension with a 
small number of points while the number of points in the 
other four dimensions goes to infinity, is legitimate. One 
could after all object that the number of lattice points in 
the standard four dimensions is also bound to be finite in 
any Monte Carlo simulation (although it can reach the 
order of 10^ for each of the coordinates, a number cer- 
tainly larger than any at present practical value of A'^^)- 

We can however observe that in the usual four- 
dimensional momentum space one does not encounter so 
wide variations for the Feynman integrand functions as 
instead is the case in the fifth dimension, where the Feyn- 
man integrands contain functions with an exponential be- 
havior, responsible for the fact that for small Ng (of the 
order of ten) large deviations of the integrals from their 
asymptotic values can be observed. Indeed, refining the 
integration grid in the usual four-dimensional momentum 
space from, say, = 60 to L^: = 80, generally produces 
minimal differences in the momentum integrals, and it is 
only in order to extract five or more significant digits that 
one needs to use = 100 or higher. I think that it is 
important in this regard to keep in mind that the shape 
of the propagator in the fifth dimension corresponds to 
physical effects strictly connected with chirality (or lack 
thereof), whereas in 4-dimensional momentum space it 
corresponds to finite size effects only. One could also take 
the point of view in which a lattice spacing for the fifth 
dimension, which is distinct from the four-dimensional 
lattice spacing a, is introduced, and imagine a situation 
where one keeps Ug finite while the usual a goes to zero. 

In the domain-wall perturbative calculations for infi- 
nite Ns made so far it was implicitly assumed that one 
could take the limit Ns — > oo for the quark propagators 
before actually using them to perform the computation of 
the Feynman diagrams and the momentum integration. 
This could in principle present some problems, because 
these asymptotic propagators know about the Atiyah- 
Singer index, which is something which is not present 
in the theory at finite Ns, and furthermore they possess 
extra infrared singularities because of the exactly mass- 
less chiral mode. In the present work we always in the 
first place compute the Feynman diagrams and the inte- 
grals using the exact Feynman rules appropriate for finite 
Ns, and only afterwards we try to investigate the limit 
Ns ^ oo hy computing the integrals anew for each in- 
creasing value of Ns- We can then observe that, apart 
from extreme values of M (see Sect. IVIIip . these inte- 
grals rapidly approach the values obtained with the the- 
ory which uses the Ns = oo propagators. By using the 
exact theory with finite Ns, we have hence also provided 
the check that in the large Ns limit one indeed recovers 
the results obtained using the simpler asymptotic the- 
ory. We can thus confirm, at least numerically, that the 
inversion of the limits was legitimate. 

In the case of the contribution of the half-circle di- 



agram of the self-energy to Eg, we have also repeated 
the whole calculation (in a general covariant gauge) 
by hand, including the evaluation of the gamma alge- 
bra and the explicit computation of the sums over the 
fifth-dimensional indices. The final expressions are very 
lengthy and it is of no interest to report them here, how- 
ever this alternative procedure provides a rather inde- 
pendent check of our calculations with respect to both 
the analytic part and the precision of the numerical in- 
tegration in the 6-dimensional code. Furthermore, since 
this time we compute the sums in the extra dimension 
analytically, this saves 2 dimensions in the numerical in- 
tegration, which is reduced to the usual 4, and this is why 
we are able to give the results for Sq with more precision 
than otherwise. Of course we cannot always employ this 
procedure, given that the other diagrams are in general 
way too complicated for a computation by hand. But this 
alternative 4-dimensional code for Sq provides another 
advantage because the dependence on Ns is now com- 
puted exactly and the computational cost is the same for 
every Ng, and thus we can think of getting results also at 
higher values of Ns, which would be too expensive in the 
6-dimensional code where the computational cost grows 
as Nl 

We conclude this Section by reminding that the tad- 
poles of course can be calculated for any Ns and AI with 
an extremely high precision, which is limited only by the 
knowledge of Zq (at present known with about 400 digits 

ill). 

VI. BILINEAR DIFFERENCES 

We now consider the calculation at finite Ns of matrix 
elements of some operators. In this Section we present 
one-loop results for the (finite) differences of chirally- 
related bilinear operators, which should become zero at 
infinite Ns, that is when chiral symmetry is fully re- 
stored. Since the vector and axial- vector currents renor- 
malize differently when chiral symmetry is broken, an 
estimate of chirality-breaking effects can indeed be given 
by how much for a given finite Ns the perturbative results 
for these currents differ from each other. The quantity 
A = Zy — Za provides such a measure of chirality break- 
ing. Moreover, one finds that Zy — Za = —{Zs — Zp)/2 
|54| . The fact that we obtain for A the same number 
with a very good precision whether we consider the vec- 
tor and axial-vector case or the scalar and pseudoscalar 
case then provides a compelling check of our calculations. 

We can successfully obtain all results for A by comput- 
ing only finite lattice diagrams, because of the following 
three facts. For one thing, the anomalous dimensions 
of the operators which are chirally related are the same 
also when one considers the subleading orders (see Sect. 

m, that is 7^-""^ = 7^-*'' and 7^^*'^ = 7^^''''. 
Furthermore, the continuum values of the finite parts, 

called i?^^ in Eq. H64|) . are also equal for these pairs of 
operators. Finally, the mismatch between lattice (at fi- 
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nite Ns) and continuum anomalous dimensions that we 
discussed at the end of Sect. IIVI is also the same for 
the chirally-related operators, and thus it cancels in their 
differences. We remind that the vector and axial- vector 
currents have now a nonvanishing anomalous dimension, 
which becomes zero only when Ng = oo. 

As we can gather in Tables IXVIIII and IXIXl the de- 
crease of the amount of chirality breaking connected to A 
follows a pattern similar to the one that we have encoun- 
tered in the case of the critical mass and to the one that 
we will see for the operator discussed in the next Section, 
that is A is rather large for small Ng and large |1 — M\, 
and decreases when Ng grows or when |1 — M| tends to- 
wards zero. Again, the part proportional to (1 — a) of 
A shows a violation of gauge invariance, since if gauge 
invariance were respected these numbers would have to 
be zero, as it happens for Wilson or overlap fermions. 
We can notice that here the deviations from gauge in- 
variance are of the same order as the ones for the critical 
mass, and that they are small. The operator studied in 
the next Section gives also very similar results regarding 
this point. 

Taking a typical case, where M = 1.8 and Ng = 16 
and ~ 2 GeV, the numbers that we have obtained 
imply for QCD a chiral violation of about 2 MeV, with 
the quenched case giving the slightly lower value. These 
numbers are much smaller, at given M and Ng, than the 
ones obtained for nireg- This could mean that the chiral 
violations are here of CMrn'L^, as it has been suggested 
for other quantities in [l7l l42|. 

It is also interesting to compare the numbers presented 
in this Section to the results that one obtains for this 
quantity when Wilson fermions are used. In this case 
one gets A = 4.82152 for the unimproved case and A = 
1.53633 for the improved case, which correspond to chiral 
violations which are at least of order 30 MeV for /3 ~ 6.0. 



VII. A POWER-DIVERGENT MIXING 

We think that it is also very instructive to study the 
case of an operator mixing which gets completely sup- 
pressed only when chiral symmetry is fully operative, so 
that the nonzero amount of residual mixing present for 
any given finite Ng provides another quantitative mea- 
sure of chiral violations. One of the simplest examples of 
this kind of mixings is probably furnished by the operator 



(72) 



which taken together with Oa2,d ~ '?7{i75^4}9 (or, in 
an alternative representation of the hypercubic group, 
Oa2,e = qiAlbD^q - 5Z]j=i97i75Ag), determines the 
first moment of the 52 structure function, which mea- 
sures the distribution of the (chiral even) transverse spin 
of quarks inside hadrons, and also receives contributions 
from twist-3 operators. More details on these operators 
can be found in [ssl Is^ (of which we follow the nota- 
tion) and references therein. Matrix elements of the op- 



erator Odi have been recently simulated using quenched 
domain- wall fermions with the DBW2 gauge action [l8l| . 
We remind that the symbol [] denotes antisymmetriza- 
tion over the relevant Lorentz indices, while {} denotes 
symmctrization, and that for the covariant derivatives 

D= D — D we use the lattice discretizations 

= 7; \u^{x)q{x + fi)-Ul{x- ii)q{x- fi)] (73) 



1 r 



q{x) D^^ - q{x 



-fi)Ul{x)-q{x-fi)U^{x-fi) . (74) 



We point out that many of the mixings which one en- 
counters in the study of physical processes, like for in- 
stance the ones occurring in the renormalization of the 
second moment of the unpolarized parton distribution 
arise as a consequence of the breaking of Lorentz 
symmetry (or, in other cases, of other symmetries apart 
from chirality) and hence do not interest us much in the 
present context. The power-divergent mixing on the lat- 
tice of Odi with an operator of lower dimension, which 
can be written as 



Cmix • — qc^iilbq-, 
a 



(75) 



is instead particularly interesting because is only caused 
by the breaking of chirality, and hence it provides a quan- 
titative measure of how much chiral symmetry has been 
broken. The finite coefficient Cmix must vanish for infi- 
nite Ng, when chiral symmetry is fully restored, and Od^ 
becomes in this case multiplicatively renormalized (like 
it also is when overlap fermions are used). For finite Ng 
the values of this coefficient in the Feynman gauge are 
reported in Table IXXI where we can see that in general 
this mixing can be considered to be almost negligible. 
Were this not the case, the removal of these lattice ar- 
tifacts in Monte Carlo simulations would become quite 
challenging. 

The results show that large chiral violations are present 
only for very small Ng, or when M is rather close to or 
2. Taking a typical case, for M = 1.8 and Ng = 16, and 
a^^ ^ 2 GeV, the chiral violations come out of about 
3 MeV, with the quenched case giving the slightly lower 
value. These numbers are much smaller, at given M and 
Ng, than the ones obtained for rrires- It could be that, as 
it has been suggested for other quantities in 0, j the 
chiral violations are here quadratic in m^es, that is this 
coefficient is doubly suppressed. 



Again, as can be deduced from Table IXXIl gauge in- 
variance is lost in the theory at finite Ng. The part pro- 
portional to 1 — a of Cmix vanishes only when Ng — 00 
(and this happens also when for example Wilson or over- 
lap fermions are used). 

Wilson fermions also suffer from this power-divergent 
mixing which is caused by the breaking of chirality, but 
the mixing coefficient is in this case gauge invariant. 
However, it takes the values Cmix = 16.243762 for the 
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unimproved and Cmix = 8.798732 for the improved case 
with Csw = 1, so that compared to domain- wall fermions 
at the standard choices of Ng and AI, the Wilson viola- 
tions are then about two orders of magnitude higher. For 
overlap fermions of course Cmix is zero. 

VIII. TOWARD THE BORDERS 

It it also interesting to see what happens when the 
domain- wall height AI is chosen dangerously close to or 
2, that is at the edges of its (at lowest order) allowed val- 
ues. After investigating in various cases, we have found 
as a general phenomenon that the exponential decay in 
Ng of the chiral violations slows down when approaching 
these borders. In some instances it can even happen that 
for a fixed M chosen very near to one of the edges the 
chiral violations increase with Ns , at least up to a certain 
point. 

One can observe this behavior already by looking at the 
tadpole diagrams, which can be calculated exactly. For 
M ^ they assume in fact the asymptotic expressions 
(in Feynman gauge) 

^MT.(iV. + 3 + A) (76) 

in the case of Sqi ^^id 

^AIT,(2N, + 3+^^ (77) 

in the case of Ei. For M — ^ 2 one just takes into account 
the fact that these tadpoles are respectively odd and even 
upon reflection around the point M = 1. We can then 
immediately see that for a given Ng the tadpoles tend 
to zero when M approaches or 2. From Table IVTl we 
can also see that the decay in Ng (at a constant AI) of 
the tadpole of Eq tends to become slower and slower as 
M nears the borders, until the rate of decay probably 
vanishes at some point and the tadpole then reaches the 
asymptotic regime of Eq. H7fci|) , where actually the chiral 
violations grow with Ng (at fixed AI). We can observe a 
similar behavior also for the other quantities considered 
in this article, which we have run for values of M close 
to |1 - AI\ - 0.99 and |1 - M| ~ 0.999. For A and c^^^ 
the rate of the exponential decay in Ng keeps decreas- 
ing when one approaches the borders, until a likely final 
disappearance, even though these quantities instead in- 
crease in M (for a fixed Ng) when M ^ and AI ^ 2. 
In fact, Cmix tends to diverge very fast when M — + 0. 

The slowing of the exponential decays of nires, A and 
Cmix when one moves AI such that |1 — M| — > 1 can be 
related to the decrease of the mobility e dge towards 
zero in these extreme regions of AI |43ll43.l45| . Choosing 
a value of M too close to or 2 can then become danger- 
ous, because the mobility edge has to remain well above 
zero in order to perform reliable Monte Carlo simulations, 
otherwise the restoration of chiral symmetry can become 



problematic. The fall of the mobility edge to zero signals 
the onset of the Aoki phase, and this can be pictorially 
seen for example in Fig. 1 of Ref. ^3 or where in 
fact it corresponds to moving, for go not too large, from 
the rightmost C phase towards the B phase, which one 
eventually enters through one of the thin "fingers" . 

We can thus briefly sum up the behavior of the chiral 
violations when M changes as follows. The coefficient 
of their exponential decay is zero or rather small when 
AI is in the vicinity of 1. It then grows when AI moves 
towards either or 2, before decreasing again and at last 
getting rather small for AI ^ and AI 2. 

A remarkable thing that we have observed for values 
of AI very near the edges is that rures and A happen 
to be smaller for small Ng than for larger values of Ng, 
that is the chiral violation initially grow in Ng instead of 
decaying exponentially. This observed initial growth is 
just a temporary one, before eventually m^es or A starts 
to decay as expected. But it could be that eventually a 
behavior like the one of Eqs. H76|l and H77|l sets in, that 
is the exponential decay disappears altogether and one 
can only see a steady increase with Ng. 

For Cmix this does not seem to happen, that is we al- 
ways observe a clear exponential-like decay, without any 
initial growth in Ng, at least for Ng > 8 and for the values 
of AI that we have investigated, up to 2 — AI = 10^^". 

Numerical examples that illustrate this phenomenon of 
the initial growth in Ng can be seen in the case of rureg in 
Tables IXXITl and IXXIIll which refer to M -> 2 (similar 
things can also be seen for the other limit, AI ^ 0). The 
onset of this behavior takes place around AI = 1.9, as 
can also be observed in the last line of Tables IXlTI to IXVl 
Furthermore, we can notice that the exponential decay 
which follows these regions of initial growth in Ng sets 
in at ever higher values of Ng when one gets nearer and 
nearer to 2 (note that for M > 1.95 the maximum of rrireg 
is reached for Ng > 28). We can then speak of a sort of 
suppression of the chiral violations for small Ng, which 
becomes stronger as AI approaches the borders, and it 
could be that in this case the density of eigenvalues or 
the radius of the modes in Eq. ^ change in such a way to 
produce this kind of effect, or that this equation breaks 
down in this region. It could also happen that the value 
of Ng for which the chiral violations reach their maximum 
(at fixed M) is further and further shifted toward higher 
values of Ng until one cannot observe any exponential 
decay at all even for very large Ng, and possibly for all 
Ng, as in the exact asymptotic results of Eqs. (|76|l and 

In Tables IXXIVI and IXXVI we can observe a similar 
behavior in the case of the difference of the vector and 
axial-vector currents. Since Cmix apparently does not in- 
stead show these effects, it could be that this quantity is 
described by a somewhat different formula than Eq. . 

It would be interesting to carry out further studies 
about these phenomena which happen near the borders 
of AI, and more investigations in the future could clarify 
these issues. We have however seen that already at one 
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loop we can observe interesting features, in some cases 
corresponding to what is expected from general consid- 
erations which can be derived from other methods. 

TABLE XXII: The residual mass near M = 2, at /3 = 5.2 in 
Feynman gauge. 



TABLE XXV: Coefficient of 5^ for A near M = 2, in Landau 



M 


Ns=8 


iV. = 12 


Ns = 16 


Ns = 20 


Ns = 24 


Ns = 28 


1.91 


0.043 


0.094 


0.120 


0.127 


0.120 


0.107 


1.92 


0.034 


0.085 


0.115 


0.128 


0.128 


0.120 


1.93 


0.026 


0.074 


0.106 


0.125 


0.132 


0.130 


1.94 


0.018 


0.061 


0.094 


0.117 


0.130 


0.134 


1.95 


0.011 


0.048 


0.080 


0.104 


0.121 


0.132 



TABLE XXIII: The residual mass near M = 2, at /? = 5.2 in 
Landau gauge. 



M 


Ns = 8 


Ns = 12 


A^. = 16 


Ns = 20 Ns 


= 24 Ns 


= 28 


1.91 


0.049 


0.099 


0.124 


0.130 


0.122 


0.108 


1.92 


0.040 


0.089 


0.119 


0.131 


0.130 


0.122 


1.93 


0.031 


0.078 


0.110 


0.128 


0.134 


0.132 


1.94 


0.023 


0.065 


0.098 


0.120 


0.132 


0.137 


1.95 


0.015 


0.052 


0.083 


0.107 


0.124 


0.134 



TABLE XXIV: Coefficient of f for A near M = 2, in Feyn- 
man gauge. 



M 


Ns = 8 


Ns = 12 


Ns = 16 


Ns =20 Ns^ 24 


Ns = 28 


1.94 


3.699 


3.151 


2.488 


1.915 


1.467 


1.125 


1.95 


4.246 


3.858 


3.220 


2.594 


2.063 


1.637 


1.96 


4.885 


4.753 


4.220 


3.592 


2.995 


2.475 


1.97 


5.630 


5.883 


5.586 


5.063 


4.474 


3.897 


1.98 


6.494 


7.301 


7.442 


7.225 


6.823 


6.332 


1.99 


7.495 


9.073 


9.945 


10.377 


10.524 


10.481 



IX. CONCLUSIONS 

In this article we have presented the calculation of a 
few one-loop amplitudes for domain-wall fermions at fi- 
nite Ns, with the intention of gathering some knowledge 
about the extent of chiral symmetry breaking for choices 
of Ns which are far away from the case of infinite ex- 
tent in the extra fifth dimension. In particular, we have 
studied three quantities whose deviations from their val- 
ues at Ns — 00 can provide some significant estimates 
of chiral violations: the residual mass, the difference be- 
tween the vector and axial-vector rcnormalization con- 
stants, and a power-divergent mixing of a deep-inelastic 
operator which is entirely due to the breaking of chirality. 

We have automated the perturbative calculations by 
developing suitable FORM codes. We have found that 
the pattern of the deviations from the case of exact chi- 
rality turns out to be approximately the same for all 



gauge. 


M 


iV. = 8 


Ns = 12 


Ns = 16 


Ns = 20 


Ns =24 


Ns = 28 


1.94 


2.112 


1.897 


1.558 


1.237 


0.973 


0.763 


1.95 


2.416 


2.310 


1.997 


1.654 


1.346 


1.090 


1.96 


2.774 


2.836 


2.602 


2.268 


1.929 


1.621 


1.97 


3.193 


3.503 


3.432 


3.179 


2.856 


2.522 


1.98 


3.682 


4.346 


4.564 


4.522 


4.335 


4.071 


1.99 


4.251 


5.402 


6.098 


6.488 


6.673 


6.718 



quantities studied, that is these violations substantially 
increase when Ns becomes small or when M approaches 
or 2, i.e., close to the borders of the region of allowed 
values of M . Our perturbative calculations show indeed 
that the numerical deviations from the case of infinite 
extension in the fifth dimension depend, apart from Ns 
and to a smaller extent from the bare coupling 1707 very 
strongly on the choice of M. These deviations can be- 
come rather pronounced when M is close to the borders 
of the region of allowed values. For M = 1.8, a stan- 
dard choice in Monte Carlo simulations, chiral violations 
remain still not small for Ns = 16, as can be seen from 
rrires, which for a lattice spacing of 2 GeV is equal to 
about 100 McV in the quenched case and about 120 MeV 
in full QCD. For the difference between the vector and 
axial-vector rcnormalization constants as well as for the 
power-divergent mixing the chiral violations are instead 
of about 2-3 MeV, suggesting that they are of higher or- 
der in rrires- 

We think that to the extent that one-loop perturbation 
theory can provide hints to the actual behavior of these 
quantities, and if higher-order corrections do not strongly 
compensate for these 1-loop effects, our work shows that 
chiral violations at Ns = 16 are still somewhat rather 
pronounced, and it suggests that one should perhaps con- 
sider performing simulations which use at least Ns = 20 
or Ns = 24 points in the fifth dimension, in order to be 
able to work with reasonably small chiral violations. 

Our work also confirms (though only from the numeri- 
cal side) the legitimacy of the assumptions made in previ- 
ous calculations at infinite Ns , where from the start it was 
postulated that the limit Ns —^ 00 for the quark propaga- 
tors could be taken before doing the actual computations 
of the Feynman diagrams and the numerical integration. 
The asymptotic propagators introduce in fact spurious 
infrared singularities not present in the theory at finite 
Ns, and in passing from Ns = 00 to any finite value of 
Ns the analytical Atiyah-Singer index of the Dirac oper- 
ator, which protects the quarks from acquiring a nonzero 
mass, disappears. One cannot then be certain that this 
transition occurs with a smooth and continuous behav- 
ior, alth oug h at the end it is likely that no problem is 
present [52|. We can observe that the results which we 
have presented in this work rapidly approach for large 
Ns the numbers which were obtained with the simplified 
asymptotic theory, with the exception of the regions very 
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close to the borders of allowed values, M — > or M ^ 2, 
where other effects can come into play. We can then a 
posteriori confirm that, at least for one-loop calculations 
made at choices of M which are not too extreme, the 
inversion of the limits was numerically legitimate. 

A disturbing finding of our calculations is the patho- 
logical behavior of renormalization factors, which are 
no longer gauge invariant, although the deviations from 
gauge invariance are not numerically large, even for val- 
ues of M far away from 1 (the pattern of chiral violations 
seems to be the same for all quantities considered in this 
work). It could be that the act of choosing a definite 
gauge affects the amount of chirality-breaking effects in 
Monte Carlo simulations as well. Furthermore, the coef- 
ficients of the logarithmic divergences turn out to depend 
on Ns and M, and only for Ng = oo they are the same as 



the ones calculated in the continuum. In particular, the 
anomalous dimension of the vector (as well as the axial- 
vector) current is not zero, and not constant in and 
M. It would be interesting in the future to investigate 
in more detail these phenomena, which could perhaps 
contribute to a fuller understanding of domain-wall sim- 
ulations. 
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